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ABSTRACT 


A  specially  designed  mass  spectrometer  was  used  to  analyze 

the  equilibrium  concentrations  of  gas  phase  solvated  ions  arising 

from:  1.  the  irradiation  of  water  vapor,  2.  the  irradiation  of 

ammonia  vapor,  and  3.  the  thermal  ionization  of  potassium  in 

water  vapor.  The  analyses  proved  the  occurrence  of  the  general 

reactions  1-3  respectively.  The  notation  n-l,n  is  used  to 

denote  specific  reactions.  (For  example,  reaction  1,2  of  the 

general  reaction  1  is  H+(IhO)  +  HO - >  H+(H„0)„). 

2  2  ■< -  2  2 

(1)  H+  (H  O)  +  H  O 5-  H+(H  0)  n  =  2-7 

2  n-1  2  ■< -  2  n 

(2)  NH  +  (NH  )  _  +  NH  - >  NH.+  (NH0)  n  =  1-6 

4  3  n-1  3  -< - 4  3  n 

(3)  K+ (HO)  +  HO  >  K+(H  0)  n  =  1-6 

2  n-1  2  < -  2  n 

Measurement  of  the  gaseous  solvent  pressure  led  to  determination 
of  the  equilibrium  constants  ^  for  reactions  1-3.  Typically 

the  pressure  was  varied  from  0.2  to  3  torr  to  prove  the  necessary 

condition  that  K  ,  is  pressure  independent.  A  temperature 

n-1  ,n 

study  of  the  reactions  over  the  range  -60°  to  600°C  allowed 

determination  of  AG°  .  ,  AH°  ,  and  AS°  from  van't  Hoff 

n-l,n  n-l,n  n-l,n 

type  plots.  With  one  exception  these  values  monotonically 
increased  with  increasing  cluster  size.  The  ranges  for  AH0 

n-l,n 

(298°K,  1  atm)  were  -36  to  -11.7,  -27  to  -7.5,  and  -17.9  to  -10.0 
kcal/mole  for  reactions  1-3  respectively.  The  experimental 
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entropy  values  for  AS0  for  reaction  1  and  As0  for 

-L  /  ^  f  -L 

reactions  2  and  3  were  in  agreement  with  those  calculated 
from  statistical  mechanics. 

The  values  of  AH0  ,  for  reaction  3  were  in  approximate 

n-l,n 

agreement  with  the  theoretical  values  obtained  from  classical 
electrostatic  calculations.  These  calculations  took  into 
account  ion-dipole,  ion-polarizability,  ion-atom  repulsion, 
and  dipole-dipole  repulsion  energy  terms.  The  calculations 
were  performed  with  the  assumption  that  all  K+  —  H^O 
distances  are  equal,  i.e.  the  cluster  has  a  single  shell  of 
water  molecules.  The  conclusion  from  these  calculations 
was  that  the  increase  in  dipole-dipole  repulsion  energy 
with  cluster  size  primarily  accounted  for  the  increase  in 

AH0  „  with  increasing  n.  A  special  electrostatic 

n-l,n 

calculation  indicated  that  the  formation  of  an  outer  shell 

of  water  molecules  becomes  energetically  favorable  for 

reaction  3  with  n  =  7  i.e.  reaction  6,7. 

The  clusters  NH  .+  (NH_)  and  H+(H  0)  are  more  complicated 

4  3  n  2  n 

since  these  clusters  could  also  be  H+(NH  )  and  HO  (HO) 

3  n+1  3  2  n-l 

respectively.  The  smooth  change  of  the  thermodynamic  values 
with  n  for  reaction  1  was  interpreted  to  mean  that  the 
appropriate  formula  for  the  hydrates  is  H+(H20)^. 

*1“ 

In  contrast  to  the  thermodynamic  values  for  H  (H20)  , 

.f 

the  results  for  NH  .  (NHJ  had  a  break  in  the  thermodynamic 

4  3  n 

values.  The  break  appeared  between  reactions  3,4  and  4,5  as 


- 


■ 
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as  a  sharp  increase  in  the  values.  This  suggested  the 

existence  of  an  inner  cluster  NH  *(NH.J Since  a  more 

4  3  4 

distantly  bound  ammonia  molecule  will  be  more  weakly  bound, 
the  values  of  AH°  and  AS0  are  greater. 

fr  ^  J  4  ^  J 

From  the  present  results,  some  comparisons  are  made 

to  ionic  solvation  in  the  liquid  phase. 

With  the  experimental  value  of  AH°  for  reaction  3 

u  f  1 

and  known  thermodynamic  values,  the  proton  affinity  of  KOH 


was  calculated  to  be  262  kcal/mole. 


* 
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1.  INTRODUCTION  TO  THE  GAS  PHASE  SOLVATION  OF  IONS 

1. 1  The  Nature  of  Solvation  Reactions 

Gas  phase  ion-molecule  reactions  have  been  classified 
according  to  reactions  1.1  -  1.5  (1): 


1.1 

Simple  Charge  Transfer 

+ 

AB 

+ 

CD 

— y 

AB  +  CD+ 

1.2 

Charge  Transfer  Induced 

Dissociation 

AB+ 

+ 

CD 

— y 

AB  +  C  ■ 

1.3 

Radical  or  Atom  Abstraction 

a. 

AB+ 

+ 

CD 

— y 

ABD+  +  C 

+ 

+ 

b. 

AB 

+ 

CD 

— y 

A  +  BCD 

1.4 

Ionic  Abstraction 

a. 

+ 

AB 

+ 

CD 

— y 

A  +  BCD+ 

+ 

+ 

b. 

AB 

+ 

CD 

— y 

ABD  +  C 

1.5 

Attachment 

+ 

AB 

+ 

CD 

— y 

ABCD+ 

Gas  phase  solvation  reactions  are  ion-molecule  reactions  of  the 
attachment  type.  More  precisely/  they  form  a  sub-class  of 
attachment  reactions.  These  reactions  are  characterized  by 
the  resulting  relatively  weak  bond  strengths  (of  the  order  of 
10  kcal/mole)  between  the  ion  and  the  molecule.  The  bonding 
between  the  ion  and  the  molecule  results  from  the  electrostatic 
attraction  of  the  ion  for  the  molecule.  In  the  present  work, 
in  which  the  solvent  molecules  were  either  NH^  or  H^O,  the 
major  attractive  forces  were  ion-dipole  and  ion-polarizability. 

Reactions  1.6  -  1.9  represent  solvation  reactions  where 
M  is  the  third  body  required  to  remove  the  excess  energy  which  is 


: 's 


' 
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due  to  the  heat  of  reaction.  If  the  cluster  has  a  large  number 


of  degrees 

of  freedom,  the 

deactivation  of 

the  cluster  by  the 

third  body 

may  not  be  necessary.  The  equilibrium  cluster  sizes 

1.6 

A+ 

M  + 

+  S  - - -*  A  S 

(0,1) 

1.7 

A+S 

M 

+  s  — ►  + 

< -  A 

2 

(1/2) 

00 

• 

,H  j 

+ 

AS2 

M  +„ 

+  S  ( - a  S3 

(2,3) 

1.9 

A+S 

n- 

_  M  + 

+  S  - - >  A  S 

1  < -  n 

(n-l,n) 

formed  during  the  solvation  reactions  are  strongly  temperature 
dependent. 

When  the  ionic  charge  becomes  well  shielded  by  one  or 
more  layers  of  solvent  molecules/  the  reactions  might  be  better 
termed  condensation  reactions,  1.10,  rather  than  solvation 


1.10 


S 


(A) 


reactions.  This  viewpoint  considers  the  cluster  A 
equivalent  to  a  liquid  drop.  In  the  present  work  it  was  found 
that  the  thermodynamic  values  for  AH° ,  AG° ,  and  AS°  for  reaction 
1.11  were  approaching  the  values  for  reaction  1.12.  The 


1. 11 

H+(H  0)  +  H_0  - *  H  (HO) 

2  n-1  2  2  n 

n 

i.  12 

H_0  .  x - *  H^0ro^ 

2  (g)  2  (£) 

comparison  to  a  liquid  drop  can  be  extended  by  considering  other 


* 
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properties  of  the  ionic  cluster.  For  example,  the  dielectric 
constant  of  the  cluster  is  expected  to  be  greater  than  that 
of  the  surrounding  gas.  As  the  cluster  size  increases,  the 
dielectric  constant  is  expected  to  increase  and  perhaps 
eventually  approach  the  dielectric  constant  of  liquid  water. 
The  dielectric  constant  of  the  cluster  is  relevant  to  the 
present  work.  If  the  dielectric  constant  increases  rapidly 
with  cluster  size,  the  strength  of  the  "electrostatic  ion- 
solvent  molecule  bonds"  will  decrease  rapidly  with  cluster 
size  because  of  the  inverse  dependence  of  the  ion-dipole 
energy  on  the  dielectric  constant. 

Ion-solvent  molecule  interactions  can  also  be  compared 
to  the  liquid  phase.  Two  terms  which  are  used  to  classify 
these  interactions  in  the  liquid  phase  are  close  and  distant 
association.  Close  association  is  the  association  due  to 
solvent  molecules  which  immediately  surround  an  ion  and  have 
their  dipoles  closely  aligned  with  the  electric  field 
gradient  of  the  ion.  Other  solvent  molecules  which  are 
farther  away  (distant  association)  are  more  free  to  rotate 
their  dipoles  away  from  the  electric  field  gradient.  Since 
the  clusters  studied  in  this  work  were  relatively  small  (up 
to  seven  solvent  molecules) ,  the  solvation  interactions  in 
these  clusters  were  probably  predominantly  due  to  close 


association. 
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1.2  Instrumentation  Used  to  Observe  Solvation  Reactions 

Mass  spectrometers  have  been  used  to  study  ion-molecule 
reactions,  including  solvation  reactions.  A  conventional  mass 
spectrometer  operates  as  follows.  In  the  mass  spectrometer 
ion  source,  an  electron  beam  is  used  to  ionize  gas  molecules  S. 

The  ions  which  are  created  include  the  parent  ion  S+  and  the  lower 
molecular  weight  ions  which  arise  from  the  fragmentation  of 
the  parent  ion.  A  positively  biased  electrode,  the  repeller, 
accelerates  the  ions  toward  the  ion  exit  slit  of  the  ion 
source.  Outside  the  ion  source  the  ions  are  further 
accelerated  and  mass  analyzed.  The  mass  analyzed  ionic 
species  are  then  individually  collected  and  recorded  as  an 

electrical  current.  When  used  in  routine  analytical  analyses, 

-5 

the  pressures  in  the  ion  source  are  in  the  range  of  10  to 

_7 

10  torr.  Under  these  conditions  the  fraction  of  ions 
undergoing  even  a  single  collision  is  negligibly  small.  Thus  in 
order  to  study  higher  order  reactions  such  as  1.6  to  1.9 
the  pressure  must  be  greatly  increased.  The  present  work 
used  a  modified  mass  spectrometer  with  ion  source  pressures 

from  0.2  to  6  torr.  Since  the  ion  acceleration,  mass  analysis, 

-4 

and  collector  regions  should  be  maintained  at  less  than  10 
torr  for  the  mass  spectrometer  to  function  properly,  the 
mass  spectrometer  used  in  the  present  work  used  a  high 
capacity  pumping  system  and  an  ion  source  with  narrow 


•- 

- 
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slits  (to  reduce  the  gas  flow  out  of  the  ion  source) . 


1.3  Thermodynamic  Equations  Pertinent  to 
Gas  Phase  Solvation  Reactions 

At  sufficiently  high  pressures  the  initial  ions  in  the 

present  work  reacted  to  form  solvated  ions.  Assuming  that 

equilibrium  is  achieved,  the  equilibrium  constant  K  ,  for 
^  n- 1 ,  n 

the  solvation  reaction  n~l,n  (reaction  1.9)  can  be  calculated 
from  equation  1.1  from  a  knowledge  of  the  partial  pressures 


1. 1 


K 


P  + 
A  S 


n-l,n 


n 


P  +  P 

A  S  ,  S 
n-1 


+  + 

of  A  S  ,  A  S  .  ,  and  S.  Since  the  partial  pressures  of 

n  n-1 

+  + 

A  S  and  AS  _  in  the  present  work  were  always  at  least 
n  n-1 

four  orders  of  magnitude  less  than  the  pressure  of  S,  to  a 
good  approximation  the  ion  source  pressure  was  equal  to  P  . 
From  the  mass  analysis  of  the  ionic  clusters,  the  ionic 
pressures  can  be  determined  from  equations  1.2  and  1.3 


where  14- 
A  S 

n 

and  1,+  are  the 

A  S 

n-1 

mass  analyzed  ion  currents 

1.2 

P  + 

=  c.  I  + 

A  S 

n 

1  AS 

n 

1.3 

P  + 

=  c_  I  +_ 

A  S  . 
n-1 

2  A  S  _ 
n-1 

+ 

of  A  S  and 
n 

A+S  ,  respectively. 
n-1 

c^  and  c^  are  proportionality 

constants  whose  evaluation  requires  both  a  knowledge  of 
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the  flow  rate  out  of  the  ion  source  and  a  knowledge  of  the 
fraction  of  each  ion  current  which  is  lost  in  transmission 
from  the  ion  source  to  the  ion  collector.  On  substitution 


of  1.2  and  1.3  into  1.1,  equation  1.4  is  obtained.  As 


suming 


1.4 


K 


c_  I  + 

1  AS 


n-l,n 


n 


c_  I  +  p 

2  A  S  S 
n-1 


°1  =  C2  reduces  equation  1.4  to  1.5.  The  validity  of  this 


1.5 


K 


IA+S 


n- 1 ,  n 


n 


I  +  P„ 

AS  _  S 
n-1 


assumption  is  further  discussed  in  Chapter  3. 

Other  thermodynamic  functions  can  be  obtained  from 

O 

K  _  values.  The  free  energy  change  AG  n  is  qiven  by 
n-l,n  ^  3  n-l,n  r  1 

equation  1.6.  Since  the  values  for  K  ,  in  the  present  work 

n-l,n 


1.6 


AG  ,  =  -  RT  in  K  , 

n-l,n  n-l,n 


were  calculated  for  a  standard  state  of  1  torr,  the  values  for 

K  ,  were  multiplied  by  760  before  substitution  into  1.6  in 

n-l,n 

order  to  obtain  the  other  thermodynamic  functions  in  a 
standard  state  of  1  atm.  A  variation  of  the  temperature 

made  it  possible  to  determine  Ah°  and  As0  .  from 

^  n-l,n  n-l,n 

equation  1.7. 


1.7  £n  K  .  =  -  AH°  AS° 

n-l,n  n-l,n  +  n-l,n 


RT 


R 
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1. 4  Significance  of  Gas  Phase  Solvation  Reactions 

A.  In  Ionic  Mobility 

The  significance  of  gas  phase  solvation  reactions 
can  be  realized  from  the  number  of  fields  in  which  these 
reactions  play  an  important  role.  One  such  field  is  ionic 
mobility  in  which  the  gas  phase  transport  of  ions  is  studied 
as  a  function  of  gas  pressure  and  electric  field  strength. 
Solvation  reactions  are  used  to  explain  the  observed  ionic 
mobilities.  For  example,  in  a  study  of  the  mobility  of 
alkali  metal  ions  in  "moist"  rare  gases,  Munson  and  Tyndall 
(  2  )  have  estimated  (as  an  upper  limit)  that  the  number  of 

-f-  -f“  **f” 

water  molecules  attached  to  Li  ,  Na  ,  and  K  are  6.2,  6.0, 
and  6.3  respectively.  The  role  of  solvation  reactions  in 
ionic  mobility  has  been  reviewed  thoroughly  by  Loeb  (  3  ) . 

B.  In  the  Liquid  Phase 

One  of  the  reasons  Munson  and  Tyndall  carried  out 
mobility  measurements  was  to  obtain  an  estimate  of  the  ion 
hydration  number  in  the  liquid  phase.  Gas  phase  solvation 
reaction  studies  can  also  be  used  to  give  an  indication  of 
the  structures  of  solvated  ions.  For  example,  it  is 
concluded  from  the  present  results  that  H+(H2°)n  is  a 
better  notation  for  the  hydrates  than  H^O  (H2°^n_i*  0ne 
advantage  of  working  in  the  gas  phase  is  that  the  solvation 
of  positive  and  negative  ions  can  be  studied  separately. 


' 
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C.  In  the  Upper  Atmosphere 

Solvated  ions  have  been  observed  in  the  upper  atmosphere. 

From  rocket  data  Narcisi  and  Bailey  (  4  )  determined  that  the 

#  *f*  -f- 

ions  H  (HO)  and  H  (HO)  account  for  35%  of  the  total  positive 

^  Z  Z 

ions  in  the  D  region  (50  -  85  Km.  altitude)  . 

D.  In  Radiation  Chemistry 

Lind  (  5  )  has  briefly  reviewed  the  role  of  ionic 
clusters  in  radiation  chemistry.  The  review  states  that  ionic 
clusters  were  postulated  by  Lind  in  the  early  1900's.  However, 
the  postulate  lost  importance  in  the  1930's  after  Eyring, 
Hirschfelder,  and  Taylor  theorized  that  free  radical  mechanisms 
could  more  satisfactorily  account  for  radiation  chemistry  results. 
More  recently,  the  review  continues,  Magee  and  Funabashi  have 
made  theoretical  calculations  which  demonstrate  the  stability 
of  ionic  clusters. 

The  increasing  amount  of  information  about  ionic  clusters 
(discussed  in  section  1.5)  may  lead  to  a  greatly  increased  use 
of  ionic  clusters  in  explaining  certain  aspects  of  radiation 
chemistry.  Anderson,  Knight,  and  Winter  (6  )  and  Bansal 
and  Freeman  (  7 )  have  recently  used  solvation  reactions  in 
some  of  their  postulated  mechanisms. 

E .  In  Gas  Discharges 

Dawson  and  Tickner  (8)  have  observed  the  ions  NH4+(NH3)^ 
in  the  ammonia  glow  discharge  while  Knewstubb  and  Tickner  (9)  have 


< 


, 
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observed  the  ions  H  (f^O)^  in  the  discharge  of  water  vapor. 


1.5  Previous  Work  on  Gas  Phase  Solvation 
Reactions  Detected  Mass  Spectrornetrically 
Because  of  the  significance  of  gas  phase  solvation 

reactions  to  the  fields  mentioned  in  section  1.4,  a  number  of 
mass  spectroscopis ts  have  delved  into  the  study  of  these 
reactions.  Yang  and  Conway  (10  )  have  investigated  reactions 

1.6  with  n  =  2-5  under  equilibrium  conditions.  The  AH° 

n-l,n 


1.6 


+ 


2n 


+  0 


2 


0 


2n+2 


values  they  obtained  ranged  from  -9.7  to  -2.1  kcal/mole.  By 

•f- 

equating  the  0  —  0  bond  energy  to  -  AH0  ,  theoretical 

4  Z  Vj  f  X 

values  of  AH0  ,  were  calculated.  These  values  were  in 
n-l,n 

satisfactory  agreement  with  the  experimental  values. 

With  the  use  of  deuterated  compounds,  Sieck,  Abramson, 
and  Futrell  (11)  have  investigated  the  mechanisms  leading  to 
the  formation  of  H+(CH_0H)o  and  H+(c„H  OH)  .  The  investigation 
used  both  a  high  pressure  mass  spectrometer  and  a  tandem  mass 
spectrometer. 

Munson  (12)  has  reported  the  high  pressure  mass  spectra 

(up  to  0.3  torr)  of  H20,  CH^H,  Cf^OCH^  and  Q^CHO.  The  ions 

H+(H  0)  ,  H+ (CH  OH)  ,  H+(CH  OCH  )  ,  and  H+(CH  CHO)  where 

2  n  3  n  33n  3  n 

__  2  for  each  species  were  found  to  constitute  a  large 
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fraction  of  the  total  ionization  of  the  respective  gases  at 
0.3  torr. 

In  a  second  paper  by  Munson  (13)  the  ions  H+(CH_NH  )  , 

3  2  n 

“I-  4“ 

H  ((CH^^NH)^  and  H  ( (CH^) ^N)  where  n  =  1-2  were  reported  to 
occur  at  high  pressures  of  the  gases  CH  NH  ,  (CH  )  NH,  and 

<J  ^  fc* 

(CH^) respectively.  In  each  gas  the  solvated  ions 
constituted  at  least  40%  of  the  total  ionization  at  a  pressure 
of  0.3  torr. 

Using  a  field  emission  ion  source  at  -195°C,  Beckey  (14) 

-f- 

was  able  to  observe  the  hydrates  H  (H^O)  with  n  as  high  as 

2  n 

ten.  By  varying  the  temperature  he  determined  the  values 
of  AH°  and  AH°  to  each  be  ebout  4  kcal/mole.  These 

-L  f  cL  Z.  /  vj 

values  are  much  lower  than  the  values  determined  in  this  work. 
It  may  be  that  the  excess  energy  acquired  by  the  ion  during 
field  ionization  led  to  non-equilibrium  results  in 
Beckey's  work. 

A  number  of  investigations  have  been  reported  from 

this  laboratory  dealing  with  solvation  reactions  occurring 

under  equilibrium  conditions.  In  references  (15,16)  the 

solvation  of  NH4+  by  NH^  is  reported.  The  competitive 

4"  4" 

solvation  of  NH^  by  NH^  and  and  of  H  by  ^0  and 

CH30H  has  been  studied  (17  ,18  respectively) .  The  hydration 
of  H+  has  also  been  investigated  (19,20).  In  reference  (21), 
the  hydration  of  the  halide  ions  F  ,  Cl  ,  Br  and  I  is 
reported. 


. 
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2.  EXPERIMENTAL 

Measurements  were  made  in  a  60° ,  15  cm  radius  mass 
spectrometer  which  had  been  specially  modified  for  the  study 
of  ion  solvation  in  the  gas  phase.  The  primary  ions  were 
generated  in  the  electrically  field  free  high  pressure  ion 
source  (up  to  5  torr)  by  a  100  kev  proton  beam  from  a  proton 
accelerator  or  by  thermal  ionization  from  a  heated  filament. 

The  ions  escaped  from  the  ion  source  by  diffusion  and  mass 
flow  into  the  high  vacuum  region  containing  the  accelerating 
electrodes.  High  vacuum  was  maintained  by  the  low  conductance 
of  gas  out  of  the  ion  source  and  by  the  high  capacity  pumping 
system. 

The  ions  were  accelerated  and  focussed  by  a  modified 
Nier  type  acceleration  system.  Detection,  after  mass  analysis, 
was  obtained  by  a  17  stage  electron  multiplier  and  a  Victoreen 
Femtometer.  The  apparatus  was  designed  and  assembled  by 
G.  Collins  (  22  ) •  Features  of  the  instrument  of  particular 
importance  to  the  gas  phase  solvation  work  are  described  in 
the  subsequent  experimental  sections. 

2.1  Ins trumental  Requirements  to  Observe  Solvation  Equilibria 

The  instrumentation  required  for  the  study  of  the  gas 
phase  solvation  of  ions  must  produce  ions,  allow  them  to  react 
at  high  pressures,  and  mass  analyze  the  reactant  and  product 
ions.  Because  of  its  ability  to  satisfy  all  three  requirements, 
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a  modified  mass  spectrometer  has  been  used  in  this  study.  A 
conventional  mass  spectrometer.  Figure  2.1,  produces  ions  by 
electron  impact  in  the  ion  source.  The  electrons  are  emitted 
by  a  hot  filament  and  accelerated  through  a  potential 
difference  between  the  ion  source  and  the  filament.  The 
electron  beam  enters  the  ion  source  through  the  electron 
beam  entrance  slit  and  ionizes  the  gas  molecules  in  its 
path  as  it  traverses  the  length  of  the  ion  source.  The 
beam  is  then  collected  in  a  Faraday  cup.  The  ions  which  are 
created  are  accelerated  out  of  the  ion  source  by  a  positive 
electric  field  between  the  repeller  and  the  ion  source. 

After  the  ions  pass  out  of  the  ion  exit  slit,  they 
are  accelerated  by  a  potential  difference,  generally  equal 
to  2000  v,  and  enter  the  magnetic  mass  analyzer  tube.  Each 
ion  is  deflected  by  the  perpendicular  magnetic  field  according 
to  its  mass  to  charge  ratio  m/e.  Each  m/e  can  be  individually 
focussed  on  the  ion  collector  by  varying  the  magnetic  field. 
This  current  is  amplified  and  displayed  on  a  recorder. 
Conventional  mass  spectrometers  are  more  fully  described  in 
reference  (23)  . 

To  investigate  gas  phase  solvation  reactions,  the  ion 
source  pressure  must  be  high.  The  required  pressure  can  be 
approximately  calculated.  To  obtain  clusters  as  large  as 
six  ligands  about  the  central  ion,  the  ion  must  undergo  six 
collisions  plus  additional  collisions  to  remove  the  excess 


- 

' 
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Figure  2.1  Partial  view  of  a  conventional  mass  spectrometer. 
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energy  produced  by  the  condensation  reaction.  The  heat  of 
solvation  per  condensation  of  one  ligand  is  in  the  range 
10-40  kcal/mole.  If  each  collision  removed  half  of  this 
excess  energy,  approximately  five  collisions  per  condensation 
would  be  necessary  to  reduce  the  energy  to  RT  for  a  temperature 
of  300 °K.  If  the  removal  of  excess  energy  is  inefficient, 
the  collision  number  will  be  higher.  For  a  distance  i  of 
0.5  cm  between  the  point  of  ionization  of  a  molecule  and  the 
ion  exit  slit,  an  ion  in  a  static  gas  undergoes  150 
collisions  at  0.1  torr.  This  is  an  approximate  calculation 
from 

„ ,  .  .  2 

2.1  Collision  No.  =  v  £ 

X  2D 

where  v  is  the  ion  velocity,  A  is  the  mean  free  path  of  the 

ion  and  D  is  the  diffusion  coefficient  (  24  ) .  The  diffusion 

2 

coefficient  was  taken  as  500  cm  /sec  as  a  representative  ionic 
diffusion  coefficient  at  0.1  torr  (25).  The  average  velocity 
was  taken  for  an  ionic  mass  of  100  at  300 °K  and  the  mean  free 
path  was  calculated  from 

2.2  A  =  . 1 _ 

.2 

tt  N  d 

where  N  is  the  gas  density  and  d  is  the  average  collision 
diameter  of  the  two  colliding  species  (  26  ) •  d  was  assumed 
to  be  5  A°  by  averaging  the  molecular  diameter  of  water, 

2.7  A°  (  27),  and  the  ionic  diameter  of  a  typical  hexahydrate. 


’ 


. 
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O 

8.0A  (28).  Since  the  ion  undergoes  approximately  150  collisions 
at  0.1  torr,  the  ionic  clusters  are  probably  in  equilibrium. 

To  test  the  equilibrium  assumption  it  is  desirable  to  be 
able  to  vary  the  pressure  by  about  a  factor  of  ten  i.e. 
to  be  able  to  increase  the  ion  source  pressure  to  1  torr. 

However,  the  mass  spectrometer  ion  acceleration  and 

-4 

mass  analysis  region  must  be  maintained  at  less  than  10  torr. 

Above  this  pressure  the  ion  beam  issuing  from  the  ion  source 

is  seriously  scattered  by  collisions.  Thus  a  pressure 

4 

differential  of  about  10  must  be  maintained  betv/een  the 
ion  source  and  the  acceleration  and  analysis  region  outside 
the  ion  source.  This  can  be  achieved  by  high  capacity  pumping 
and  by  reducing  the  leakage  of  gas  out  of  the  ion  source. 

2 . 2  The  High  Capacity  Pumping  System 

In  the  present  instrument  the  total  area  of  the 
ion  source  slits  has  been  reduced  by  about  a  factor  of  100  and 
the  pumping  speed  at  the  ion  source  increased  by  about  100 
over  a  conventional  mass  spectrometer.  The  gas  was  pumped 
through  a  six  inch  liquid  nitrogen  baffle  and  a  six  inch 
gate  valve  by  a  1400  L/sec  diffusion  pump.  The  manufacturer 
rates  the  pumping  speed  at  the  entrance  of  the  gate  valve  as 
450  L/sec  (  29)  for  air.  -The  reduction  in  pumping  speed  is 
due  to  the  limited  conductance  of  the  baffle  and  valve. 
Conductances  in  series  add  according  to  equation  2.3 
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2.3 


1 

F 


l  L_ 

i  F . 

l 


where  F^  is  the  total  conductance  and  F.  is  the  conductance 
t  l 


of  each  component.  The  reduced  pumping  speed  S  is  given  by 
2.4  S 


Ft  x  Sp 


Ft+  Sp 

where  Sp  is  the  speed  at  the  diffusion  pump  (  30  )  .  By  com¬ 
parison  to  the  manufacturer ' s  rated  speed,  the  pumping  speed 
at  the  ion  source  was  measured  as  350  L/sec  for  ammonia  with 
the  liquid  nitrogen  baffle  uncooled.  The  speed  is  calculated 
from  the  pressure  drop  in  the  gas  handling  plant  reservoir 
with  respect  to  time  and  from  the  pressure  in  the  vacuum 

housing  (measured  by  an  ion  gauge) .  In  this  measurement  the 

-4 

ion  gauge  pressure  was  2.1x10  at  an  ion  source  pressure  of 

3  torr.  Cooling  the  baffle  with  liquid  nitrogen  decreased 

-5 

the  ion  gauge  pressure  to  3.8x10  torr.  This  indicates  a 
pumping  speed  increase  to  1900  L/sec  due  to  condensation  of 
ammonia  on  the  baffle. 

From  a  knowledge  of  the  pumping  speed  and  the  required 
pressure  differential  between  the  ion  source  and  high  vacuum 
region  outside  the  ion  source,  the  maximum  gas  conductance 
of  the  ion  source  can  be  calculated.  From  the  equation 


2.5 


PXF  =  P2S 


where  P  is  the  ion  source  pressure,  F  is  the  conductance  of 
the  ion  source  exit  slits,  P2  is  the  vacuum  housing  pressure. 
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and  S  is  the  pumping  speed  in  the  vacuum  housing  ( 30  )  ,  the 

3 

maximum  conductance  is  190  cm  /sec.  By  equation  2.6 

2.6  F  =  v  A 

4 

where  v  is  the  average  velocity  and  A  is  the  area  of  the  slits 

-2  2 

(  30  )  ,  the  maximum  total  area  of  the  slits  is  about  10  cm  . 
The  total  area  of  the  slits  used  was  about  an  order  of 
magnitude  smaller.  Because  the  ion  current  is  also 
proportional  to  the  slit  area,  the  slit  area  cannot  be  too 
small  lest  no  ion  current  be  detected  at  the  ion  collector. 

2 . 3  The  Proton  Generator 

A.  Principles  of  Operation 

Another  difference  in  the  present  instrument 
compared  to  a  conventional  mass  spectrometer  was  that  a  100 
kev  proton  beam  was  used  to  ionize  the  gas  in  the  ion  source. 
The  present  mass  spectrometer  was  equipped  with  a  proton  beam 
generator  for  two  reasons.  First,  the  ionization  of  atoms 
and  molecules  by  fast  protons  could  be  studied.  This  study 
was  completed  by  Collins  (  22  ) .  Secondly,  a  high  energy 
beam  is  required  to  penetrate  through  the  high  pressures  of 
gas  in  the  ion  source.  At  one  torr  ion  source  pressure,  a 
100  kev  proton  loses  3  kev  through  collisions  in  the  ion 
source  (  31  ) .  These  collisions  result  in  ionizations  and 
scattering  of  the  proton  beam.  Use  of  a  lower  energy  beam 
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would  lead  to  greater  scattering.  This  means  that  some  of 
the  ionization  would  occur  closer  to  the  ion  exit  slit  and 

that  some  of  the  ions  might  not  reach  chemical  equilibrium 
before  passing  out  the  ion  exit  slit. 

A  Cockcroft-Walton  type  proton  generator,  model 
number  100-lH,  Texas  Nuclear  Corp. ,  was  used  to  supply  a 
beam  of  100  kev  protons.  Because  the  unit  is  described 
elsewhere  (  32  ) ,  only  a  brief  description  is  given  here. 
Figure  2.2  shows  the  proton  generator  and  mass  spectrometer. 
The  protons  are  created  in  a  hydrogen  discharge  in  a  pyrex 
ion  bottle  shown  in  Figure  2.3.  The  pressure  in  the  ion 
bottle,  10-100  millitorr,  is  regulated  by  a  variable  leak. 

A  radiofrequency  field  applied  to  the  two  excitor  rings 
ionizes  the  hydrogen.  An  extraction  electrode,  o-5  kv 
positive  with  respect  to  the  accelerating  potential  repels 
the  positive  ions  out  of  the  bottle  and  attracts  electrons. 
The  electrons  cause  ionization  and  sustain  the  discharge. 

A  longitudinal  magnetic  field  supplied  by  a  solenoid  causes 
the  electrons  to  spiral  and  ionize  more  molecules  than  the 
straight  path. 

Upon  leaving  the  ion  bottle  the  protons  are  focussed 
by  a  focus  electrode,  0-10  kev  negative  to  the  accelerating 
potential.  The  beam  travels  through  the  accelerator  (a 
series  of  ten  electrodes) ,  a  magnetic  field  from  a 
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to  gate  valve,  6"  diffusion  pump. 
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perinanent  ring  magnet,  a  meter  long  field  free  drift  tube,  and  two 
collimation  slits  before  reaching  the  mass  spectrometer  ion  source 
The  magnetic  field  of  the  ring  magnet  which  is  perpendicular 
to  the  proton  beam  slightly  deflects  the  beam.  Adjustment 
of  the  ring  magnet  is  used  to  position  the  beam  onto  the 
ion  source  proton  beam  entrance  slit.  The  purpose  of  the 
collimation  slits  is  to  reduce  the  beam  intensity  incident 
on  the  ion  source.  The  heat  which  must  be  dissipated  from 
the  beam  is  the  product  of  the  proton  current  times  the 
proton  accelerating  voltage.  Since  the  manufacturer  claims 
a  maximum  current  of  1500  yA  ( 33  ) ,  a  maximum  150  W  has 
to  be  dissipated.  Since  heat  loss  is  by  radiation  and 
conduction  only,  a  collimator  could  heat  up  considerably. 

The  first  slit  was  water  cooled.  The  slit  was  located 
about  20  cm  from  the  ion  source  and  had  dimensions  of 
3mmx5mm.  The  potential  of  the  slit  was  earth  potential. 

The  second  slit,  which  was  positioned  about  2  cm  from 
the  ion  source,  had  dimensions  of  0.5mmx3mm.  Its  potential 
was  +225  v  to  the  ion  source  potential.  This  potential 
collected  the  secondary  electrons  which  were  emitted 
when  the  proton  beam  impinged  on  the  ion  source.  The 
presence  of  secondary  electrons  in  the  ion  source  makes 
the  ionization  of  the  gas  in  the  source  more  conplex. 

After  passing  through  the  ion  source  the  beam  is  collected 
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in  a  Faraday  cup.  Ion  source  geometry  will  be  discussed  in 
a  later  section. 

B.  Operating  Procedure 

The  effect  of  the  proton  beam  controls  on  the  proton 
beam  current  collected  by  the  Faraday  cup  in  the  ion  source 
is  shown  in  Table  2.1.  The  extraction  voltage  could  be 
varied  over  the  range  0-5  kv  while  the  focus  potential 
could  be  adjusted  from  0  to  -10  kv.  The  potentials  are 
relative  to  the  proton  accelerating  potential  of  0  -  100  kv. 
Although  the  accelerating  potential  was  variable  from 
0  -  100  kv,  all  the  present  work  used  100  kev  protons. 

The  data  in  Table  2.1  indicates  the  extraction  potential 
was  most  important  in  obtaining  high  proton  beam  currents  at 
the  ion  source.  The  focus  and  solenoid  controls  could  be 

left  at  zero  volts  and  at  zero  current  respectively.  Changing 

-5 

the  pressure  over  the  range  1  -  2  x  10  torr  as  measured 
by  the  Phillips  cold  cathode  gauge  varied  the  proton  beam 
current  by  a  factor  of  two. 

The  stability  of  the  proton  beam  as  indicated  by  the 
change  in  mass  analyzed  ion  current  with  time  was  normally 
about  ±5%  over  a  10  minute  period.  Generally  about  10 
minutes  was  required  to  scan  over  and  measure  the  required 
mass  analyzed  ion  currents.  Because  the  stability  was 
sometimes  poorer,  the  first  ion  current  measured  was 
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TABLE  2 . 1 


Effect  of  Proton  Beam  Controls  on  the  Proton 


Current  Striking  the  Faraday  Cup  of  the  Ion  Source 


I.  Maximum  Beam  Conditions 


Control  ‘ 
Pressure 


Extraction 


Control  Setting  Faraday  Cup  Current 

-5 

1.2x10  torr  (on  the 
downstream  Phillips 
cold  cathode  gauge) 

+  3900  v 


Focus 


800  v 


Solenoid 


0  A 


1.9xl0"?  A 


II.  Effect  of  Control  Settings  at  1.2x10 


-5 


torr  Pressure 


Focus 

— 

3200 

V 

Solenoid 

0 

A 

Extraction 

+ 

3400 

V 

0. 

,  2x10 

+ 

3600 

0. 

,7 

+ 

3800 

0. 

,  8 

+ 

4000 

0. 

.8 

+ 

4200 

0. 

,4 

Extraction 

+ 

3800 

V 

Solenoid 

0 

A 

Focus 

0 

V 

1. 

,  7 

- 

400 

V 

1. 

,7 

- 

800 

1. 

,9 

- 

1200 

1. 

,6 

— 

1600 

1. 

,  4 
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TABLE  2.1  (continued) 


Control  Setting  Faraday  Cup  Current 


Extraction 

+  3800  v 

Focus 

-  800  v 

Solenoid 

0.0  A 

1.9 

0.3 

1.9 

0.6 

0.2 
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always  remeasured  after  the  other  ion  currents  were  measured. 

If  the  two  currents  did  not  agree  within  ±5%,  the  mass 
spectrum  was  remeasured.  Any  change  in  the  ion  source 
such  as  temperature  or  pressure  would  also  have  affected 
the  ion  current. 

C.  Proton  Beam  Deflection  Plates 

To  permit  measurement  of  the  ion  source  residence  time 
of  an  ion,  deflection  plates  were  installed  in  the  proton 
beam  drift  tube.  A  pulsed  deflection  of  the  proton  beam 
coupled  with  a  pulsed  detection  of  the  mass  analyzed  ions 
by  a  variable  time  delay  can  be  used  to  determine  the  ion 
source  residence  time  of  an  ion.  From  the  initial  change  of 
ion  current  with  time,  the  forward  reaction  rate  constant 
can  be  determined.  When  equilibrium  is  reached,  the 
relative  ion  currents  become  invariant  with  time.  This 
allows  calculation  of  both  the  equilibrium  constant  and  the 
reverse  reaction  rate  constant.  The  pulsing  technique  has 
been  applied  by  many  researchers  (  33  )  at  lav  pressures  to 
determine  forward  rate  constants.  At  high  pressures  the 
technique  has  been  applied  to  gas  phase  solvation  reactions 
by  Durden,  Good  and  Kebarle  (  34  ) . 

The  design  of  the  deflection  plates  was  based  on  equation 
2.7  (the  equation  can  be  readily  derived  from  fundamental  laws  and 

2‘7  y  =  VDef 1.  *  X 

2V.  Y 

ion 

equation  11.5  of  Ref.  35). 
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where  y  is  the  deflected  distance  at  a  target,  X  is  the 

distance  from  the  plates  to  the  target,  Y  is  the  distance 

between  the  plates,  l  is  the  length  of  the  plates,  V. 

ion 

is  the  voltage  used  to  accelerate  the  ion  to  the  deflection 

plates  and  __  is  the  potential  difference  between  the 
^  Defl. 

plates.  The  plates  were  made  1.2  cm  apart  and  19.6  cm 

long.  The  center  of  the  plates  was  40  cm  from  the  ion  source. 

The  calculated  deflection  which  was  3x10  ^cm/V  agreed 

Defl. 

with  the  estimated  experimental  deflection.  To  observe  the 
deflection,  the  second  collimation  slit,  which  is  located 
2  cm  in  front  of  the  ion  source,  was  coated  with  a  phosphor. 

The  proton  beam  image  on  the  collimation  slit  was  observed 
through  a  plexiglass  window. 

To  make  the  voltage  required  to  deflect  the  beam  off 
the  ion  source  entrance  slit  small  enough  to  use  a 
conventional  pulsing  unit,  the  first  collimation  slit  was 
moved  ahead  of  the  deflection  plates.  The  dimensions  of 
the  slit  were  reduced  to  0.25  mm  x  10  mm. 

The  observed  image  of  the  first  collimation  slit  on 
the  second  collimation  slit  had  a  width  of  about  0.5  mm. 

Applying  a  static  deflection  voltage  of  50v  not  only 

deflected,  the  beam  but  also  approximately  doubled  the 

slit  image.  In  spite  of  this  anomalous  spread  in  the  width  of  the 

proton  beam,  a  50  v  pulsed  deflection  should  be  sufficient 

to  produce  a  pulsed  beam  in  the  ion  source. 


? 
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D.  X-ray  Hazard 

If  secondary  electrons  ejected  by  the  proton  beam  at 
the  first  collimation  slit  are  incompletely  suppressed,  they 
may  enter  the  stray  positive  field  of  the  accelerator  and  be 
accelerated  up  to  100  kev.  Vflien  these  fast  electrons  strike 
a  surface.  X-rays  are  ejected.  This  is  a  normal  hazard  of 
the  accelerator.  Since  the  point  of  emission  of  secondary 
electrons  was  moved  much  closer  to  the  accelerator  with  the 
installation  of  the  deflection  plates,  a  radiation  monitor 
was  used  to  remeasure  the  X-ray  flux.  The  readings  were 
similar  to  those  taken  previously.  The  meter  read  less 
than  0.01  mrh  with  the  monitor  facing  the  lead  shield  which 
surrounds  the  accelerator.  The  lead  shield  was  an  1/8 
inch  thick  and  6  feet  high.  Directly  above  the  shield  the 
monitor  read  0.5  mrh.  The  weekly  tolerable  dose  is  100 
mr  ( 36  )  .  As  long  as  the  accelerator  is  functioning 
properly,  X-rays  should  not  be  a  hazard. 

2 . 4  Source  of  Potassium  Ions 
A.  Introduction 

Although  the  proton  beam  conveniently  produced  ions 
for  the  study  of  H+(Ho0)  and  NH.+(NH  )  ,  its  use  for  the 
generation  of  K+  proved  to  be  impractical  due  to  the  reaction 
between  water  vapor  and  potassium  vapor.  For  this  reason 
it  was  necessary  to  generate  potassium  ions  in  a  source 
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external  to  the  ion  source.  An  alternative  method  to  ionization 
by  a  high  energy  beam  is  thermal  ionization.  The  method 
depends  on  the  fact  that  evaporation  from  a  surface  may  occur 
as  the  evaporation  of  a  neutral  species  or  a  positive  ion. 

The  relative  probability  of  evaporation  as  a  positive  ion 
is  given  by  the  Langmuir- Saha  equation  (  37) 

2 . 8  n_^  rc  exp  ({>-1 

n  kT 

o 

where  n  is  the  number  of  positive  ions,  n  is  the  number  of 
+  o 

neutral  particles,  $  is  work  function  of  the  surface,  and 
I  is  the  ionization  potential  of  the  neutral  particle.  Because 
(J)-I  >>  kT,  thermal  ionization  sources  operate  at  high 
temperatures.  The  sources  which  are  used  in  mass  spectrometric 
inorganic  analyses  operate  at  2000°K  (23  ).  The  difficulty 
in  using  a  high  temperature  source  of  ions  is  that  the 
source  has  to  be  sufficiently  distant  from  the  ion 
source  so  that  the  ion  source  does  not  heat  up.  (It  was 
desired  to  study  solvation  equilibrium  at  temperatures  as 
low  as  250 °K) . 

From  a  literature  search  for  a  suitable  thermal  ionization 
source  and  from  some  experimentation,  it  was  decided  to  use  a 
Kunsman  source  located  in  a  housing  which  was  gas  tight  to 
the  ion  source.  See  Fig.  2.4.  The  thermal  ionization  in 
the  Kunsman  source  arises  from  a  specially  prepared  filament. 

This  set-up  had  several  advantages.  Because  the  operating 
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temperature  was  only  about  800°K/  the  filament  could  be 
placed  very  close  to  the  ion  source.  This  meant  that  the 
electric  field  required  to  transport  a  sufficient  supply  of 
ions  to  the  ion  source  could  be  low.  Consequently  the 
number  of  collisions  in  the  ion  source  required  to  reduce 
the  energy  of  the  ion  to  thermal  equilibrium  is  less.  By 
comparison,  the  alternative  to  this  arrangement  is  the  placement 
of  the  Kunsman  source  at  a  distance  from  the  ion  source.  The 
ions  would  then  have  to  be  accelerated  through  vacuum  into 
the  ion  source  through  an  entrance  slit.  Although  this  method 
reduces  the  heating  of  the  ion  source,  experimentally  this 
method  was  found  to  deliver  an  insufficient  ion  current  to  the 
ion  source.  The  reason  is  that  the  entrance  slit  had  to  be 
very  small  to  maintain  the  necessary  high  vacuum  in  the 
vacuum  housing  as  discussed  in  Section  2.1.  Even  if  a 
sufficient  ion  current  entered  the  ion  source,  the  potential 
energy  of  the  incoming  ion  would  have  been  an  order  of 
magnitude  higher  than  with  the  present  arrangement. 

B.  Manufacture  of  a  Kunsman  Thermionic  Source 

The  manufacture  of  a  Kunsman  thermionic  source  is 

given  in  a  review  by  Jenkins  and  Trodden  (  37).  The  first 

step  was  the  fusion  of  a  mixture  consisting  of  98%  Fe^^, 

1%  Al  0  ,  and  1%  KNO  .  The  KNO  can  be  replaced  by  any  of 
2  3  8  -j 

the  alkali  nitrates  in  order  to  generate  other  alkali  ions. 
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In  the  fused  product  the  Fe^O„  was  oxidized  to  Fe  0  .  The 

2  3  3  4 

product  was  ground  to  a  powder  before  it  was  fused  onto  a 
Pt  filament.  The  filament  was  made  from  a  strip  of  Pt  1  mm 
wide  and  2  cm  long  which  was  cut  from  a  0.001  inch  thick 
Pt  foil.  The  width  should  be  uniform  or  the  filament  will 
have  hot  spots  in  it  when  the  heating  current  is  passed 
through.  The  Pt  strip  was  spot  welded  to  two  nichrome 
studs  which  were  spaced  about  1  cm  apart.  Fusion  of  the 
powder  onto  the  filament  required  special  care  because  the 
melting  point  of  the  powder  was  close  to  the  melting  point 
of  Pt.  A  layer  of  powder  1  mm  high  was  first  spread 
uniformly  along  the  length  of  the  filament.  A  torch 
provided  with  a  special  3/4  mm  hole  tip  was  adjusted  to 
give  a  fine  intense  flame.  The  powder  was  heated  until  it 
just  began  to  melt  and  flow.  Cycles  of  melting  and  cooling 
were  used  to  prevent  the  Pt  from  melting.  If  the  coating 
is  not  uniform  over  the  filament  surface,  the  surface 
temperature  will  be  uneven  and  the  ion  emission  reduced. 

Reduction  of  the  surface  was  the  next  step.  The 
filament  was  placed  in  a  glass  chamber.  The  chamber  was 
evacuated  and  filled  with  3-10  torr  H  .  The  filament  wattage 
was  increased  until  the  light  emitted  by  the  filament  was 
just  visible.  The  necessary  wattage  was  15-30W  depending 
on  the  pressure.  The  surface  was  sufficiently  reduced  in 
about  30  min.  A  longer  reduction  period  did  not  improve  the 


filament  emission. 
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After  the  reduction  the  filament  wattage  was  decreased 
and  the  glass  chamber  pumped  out.  Increasing  the  filament 
wattage  brought  an  onset  of  positive  current  at  about  10W. 

The  emitted  current  was  collected  on  a  plate  with  a  negative 
bias.  Increasing  the  wattage  until  the  filament  became 
dull  red  yielded  a  positive  current  of  1-2  yA.  If  the 
filament  temperature  was  increased  to  a  bright  red,  a 
maximum  of  10  yA  could  be  obtained.  This  high  current 
however,  lasted  only  about  20  min.  after  which  the  current 
dropped  well  below  1  yA.  A  high  temperature  during  reduction 
was  also  found  to  be  detrimental. 

C.  Utilization  of  the  -Kunsman  Source 
Figure  2.4  illustrates  Sources  I  and  II,  the  Kunsman 
sources  used  in  the  present  work.  The  two  sources  were 
necessary  to  prove  that  the  reactions  were  not  in  thermal 
disequilibrium  due  to  the  heating  of  the  ion  source  or  the 
gas  in  the  ion  source  by  the  filament.  (The  effect  of  the 
different  sources  on  the  results  is  discussed  in  Section  5.4) 
Source  II  was  similar  to  Source  I  except  for  location. 

The  heater  current  to  the  filament  was  controlled  by 
two  variable  autotransformers  in  series  and  regulated  by  a 
constant  voltage  transformer.  The  constant  voltage  transformer 
received  110  V.A.C.  from  a  3000  v  isolation  transformer  which 
was  plugged  into  the  mains,  110  V.A.C.  The  power  required 


. 
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1  cm 

Source  I  with  5mm  orifice 


Figure  2.4  Kunsman  Sources  I  and  II  of  K  shown 


installed  in  the  ion  source. 
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to  give  a  sufficient  emission  of  ions  from  the  filament  at 

water  vapor  pressures  up  to  1  torr  was  about  20W. 

Because  the  ions  which  were  emitted  by  the  heated 

filament  diffuse  rapidly  to  the  wall,  it  was  necessary  to 

attract  the  ions  to  the  ion  source  by  a  potential 

difference  to  obtain  a  detectable  ion  current  after  mass 

analysis.  The  potential  difference  was  increased  until  the 

ion  current  after  mass  analysis  was  just  adequate  for 

accurate  measurement.  A  potential  difference  of  10  v  was 

found  to  give  a  sufficient  ion  current  at  all  pressures 

in  the  ion  source.  Subsequently  this  potential  was 

fixed  at  10  v  or  about  lOOv/cm  for  Sources  I  and  II. 

After  the  K+(H  0)  study  was  completed,  research 

2  n 

“f*  ■if- 

was  begun  on  the  K  (NH_)  system.  Unlike  the  K  (ELO) 

3  n  2  n  . 

•f* 

results,  the  K  results  were  affected  by  a  variation 

in  the  filament  wattage.  To  investigate  the  K  {NH  ) 

3  n 

•f* 

system.  Source  III  was  made.  The  K  (thO)  results  were 

2  n 

then  rechecked  at  30°C  and  found  to  be  in  agreement  with 

t 

the  previous  results.  Figure  2.5  shows  the  design  of 
Source  III.  The  filament  v/as  moved  from  1  mm  from  the 
ion  source  to  20  mm.  A  ceramic  insulator  thermally  and 
electrically  isolated  the  filament  housing  from  the  ion 
source.  Heat  energy  from  the  filament  was  reflected  by  a 
nichrome  reflector  backed  by  a  copper  plate  which  dissipated 
heat  energy  to  the  ion  source. 


« 

. 


1  cm 


Figure  2.5  Kunsman  Source  III  of  K  shown  with 
modified  Ion  Source  Top. 
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Because  the  filament  was  farther  from  the  ion  source, 
the  problem  of  obtaining  a  sufficient  steady  state  concentration 
of  potassium  ions  in  the  ion  source  was  more  difficult. 

Because  the  ion  source  was  farther  from  the  filament  than  the 
filament  housing,  the  potential  of  the  ion  source  had  to  be 
sufficiently  lower  than  the  filament  housing  potential  to 
preferentially  attract  the  potassium  ions.  To  repel  the 
ions  diffusing  away  from  the  ion  source,  a  repeller  electrode 
was  installed.  Experimentally,  the  highest  K+  currents 
to  the  ion  source  occurred  with  the  repeller  electrode, 

.  filament  housing,  and  filament  at  the  same  potential.  A 
repulsive  potential  on  the  repeller  electrode  evidently 
lowered  the  emission  of  the  filament  while  an  attractive 
potential  on  the  filament  housing  collected  ions  in  the 
long  narrow  channel  to  the  ion  source.  A  potential  of 
225  v  was  necessary  to  attract  a  sufficient  number  of  ions 
per  sec.  into  the  ion  source.  The  electric  field  strength 
was  approximately  the  same  as  the  electric  field  strength 
as  used  with  Sources  I  and  II. 

2. 5  The  Ion  Source 

The  ion  source  in  which  the  potassium  ion  clustering 
reactions  occurred  is  shown  in  Figure  2.4.  The  purpose  of 
the  ion  source  was  to  provide  a  high  pressure  ion  reaction 
chamber  at  a  uniform  temperature  which  could  be  varied  over 
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a  wide  temperature  range.  The  high  pressure  requirement  was 

discussed  in  section  2.1.  In  that  section  the  maximum  total 

-2  2 

area  of  the  ion  source  slits  was  calculated  to  be  10  cm  . 
Because  the  conductance  of  a  slit  increases  when  the  mean 
free  path  of  the  gas  is  smaller  than  the  slit  width,  it  was 
desired  to  make  the  slit  width  less  than  30  y,  the  mean  free 
path  of  water  at  1  torr,  273°K.  (  38  )  (The  increase  in 
conductance  is  due  to  a  transition  from  molecular  flow  to 
viscous  flow  (  39  ) . 

One  difficulty  in  making  such  narrow  slits  is  that  the 
thickness  of  foil  at  the  slit  edge  must  be  less  than  the  slit 
width  or  the  ions  will  discharge  themselves  on  the  wall  of 
the  slit  in  passage  through  the  slit.  As  a  result  the  knife 
edges  of  razor  blades  were  used  to  make  the  slits.  The 
knife  edges,  after  having  been  cut  to  the  appropriate  size, 
were  placed  over  a  slot,  1  mm  x  4  mm,  which  was  in  the  center 
of  a  support  which  bolted  into  the  ion  source. 

To  make  the  slit,  one  half  of  the  slit  was  spot  welded 
on  the  center  line  of  the  plate.  The  second  half  was  placed 
in  approximately  the  correct  position  and  clamped  in  the 
spot  welder.  With  a  strong  light  behind  the  slit,  the 
second  half  was  adjusted  until  the  two  knife  edges  were 
parallel  and  the  light  coming  through  the  slit  appeared 
as  a  tiny  crack. 

When  the  geometry  seemed  satisfactory,  a  single  spot 
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weld  was  made.  If  the  geometry  was  satisfactory  as  viewed 
under  a  microscope  which  could  measure  distance  to  ± 1 y ,  the 
spot  welding  was  completed  and  the  slit  rechecked  under  the 
microscope.  The  greatest  difficulty  in  the  slit  construction 
was  vertical  displacement.  Often  the  slit  had  to  be  rejected 
because  the  plane  of  one  slit  was  as  much  as  30y  above  the 
other.  The  microscope  could  also  measure  depth.  The 
dimensions  of  the  ion  exit  slit  were  generally  about  3y  x  4mm. 

The  finished  ion  exit  slit  was  then  bolted  into  the 
ion  source.  In  the  case  of  the  experiments  utilizing  the  proton 
beam,  the  ion  source  also  had  a  proton  beam  entrance  slit  and 
collimation  slit  as  illustrated  in  Figure  2.6.  These  slits 
were  first  made  from  Ta  foil.  However,  after  it  was  found 
that  the  Ta  foil  reacted  with  water  vapor  at  high  ion  source 
temperatures,  the  slits  were  made  from  razor  blades.  The 
inner  slit  was  0.4  mm  x  5  mm  and  the  outer  slit  was  0.015  mm  x 
5  mm.  The  total  area  of  the  ion  source  slits  was  thus  an 
order  of  magnitude  less  than  the  maximum  permissible 

calculated  in  section  2.1.  Because  the  vacuum  housing  pressure 

-5 

as  measured  by  an  ion  gauge  was  about  3x10  torr  at  1  torr 
ion  source  pressures  at  298°K,  the  ion  source  could  be 
operated  at  this  high  pressure  without  serious  scattering 
of  the  ion  beam  issuing  from  the  ion  exit  slit. 

Besides  operation  at  high  pressure,  the  ion  source 
had  to  have  a  uniform  temperature.  To  obtain  a  uniform 
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Figure  2.6a  Ion  Source  used  with  proton  beam.  Side  view. 


Proton  Beam 


a> 

r — I 

o 

K 

p 

0) 

ft 

fd 

<u 


a) 

i — I 

ft 

o 

o 

o 

o 

E 

p 

a) 

.c 

Eh 


> 

O 

in 

co 


> 

o 

o 

in 

ft 


ft 

P 

o 

g 

P 

Q 


> 

o 

CO 

I — I 

I 

o 

CD 


> 

o 

o 


+ 

co 

<D 

ft 

fd 

ft 

ft 

ft 

i — I 

fd 

A 


w 

p 

o 

o 

ft 


a;  <u 
o  ft 
p  nj 


P 

o 

co 


ft 

ft 


ft 
p 

o 
ft 
ft 
p  p 
O  CO 


c 

o 

H 


> 

o 

o 


> 

o 


co 
<u 
-P 
id 
■ — i 

> 

O  4h 

•  i — I 

o  id 

.c 

to 

p 

ft  p 
g  0 

ft  ft 
ft  c 
<d  <d 
P  u 

e  6 

?d  td 

a)  a; 
M  « 

4-  + 


> 

o 

• 

o 

p 

0) 

N 

ft 

ft 

id 

g 

id 

cn 

CO 

id 

e 

o 

ft 

a> 

o 

P 

fd 

p 

ft 

g 

w 


CO 

<d 

ft 

fd 

ft 

ft 

g 

o 

ft 

ft 

rd 

p 

a) 

i — I 

CD 

O 

O 

C 

g 

id 

a) 

o 

p 

P 

o 

CO 

c 

o 

H 


& 

VD 

• 

CN 

CD 

P 

P 

ft 

•H 

ft 


id 

CD 

CN 


Ion  Source  Side  View, 90°  rotated  from  side  view  in  Figure 
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temperature  the  ion  source  bottom  was  machined  from  a  single 
stainless  steel  block.  The  heater  holes  were  placed  1.5  cm 
from  the  ion  source  cavity  to  allow  the  heat  flow  to  the 
cavity  wall  to  be  uniform.  The  ion  source  tops  were 
carefully  machined  to  give  a  good  fit  to  the  ion  source 
bottom.  The  contact  between  the  ion  source  top  and  bottom 
was  improved  by  securely  fastening  the  two  parts  together 
with  four  heavy  bolts.  The  temperature  control  and  measurement 
is  discussed  in  sections  2.9  and  2.10. 

2 . 6  Ion  Acceleration ,  Mass  Analysis  and  Measurement 

of  Ion  Current 

In  the  high  pressure  ion  source,  ions  arising  either 
from  the  proton  beam  or  the  Kunsman  source  reacted  to  form 
ionic  clusters.  After  these  ionic  clusters  passed  out  of 
the  ion  source  through  the  ion  exit  slit,  they  were 
accelerated  from  the  ion  source  potential,  +1850  v,  to 
earth  potential.  Figure  2.6  shows  the  electrode  system, 
which  was  a  modified  Nier  design.  The  first  electrode 
spacing  to  the  ion  source  was  increased  from  1  mm  to 
10  mm.  As  discussed  in  detail  in  section  3.4,  the  shorter 
distance  prevented  the  space  around  ion  exit  slit  from 
being  adequately  pumped.  The  resulting  high  pressure 
caused  the  ioni'c  clusters  to  undergo  disequilibrium 


collisions  outside  the  ion  source. 
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After  acceleration  the  ions  entered  the  magnetic 
analysis  tube  and  were  deflected  60°  to  the  ion  collector. 
The  magnetic  analyzer  tube  entrance  slit  and  exit  slit 
were  0.500  inch  x  0.010  inch  and  0.500  inch  x  0.020  inch 
respectively.  The  theoretical  resolution  is  given  by 

2 . 9  =  _ R _ 

A  m  d,+d„ 

1  2 

where  m  is  the  highest  mass  which  is  fully  resolved  from 
m+1,  R  is  the  radius  of  curvature  of  the  analyzer  tube, 

15  cm  in  the  present  instrument,  and  d^  and  d^  are  the 
analyzer  entrance  and  exit  slit  widths  (40  )  .  While  the 
calculated  resolution  was  200  for  the  present  instrument, 
the  observed  resolution  was  just  slightly  less. 

After  the  ions  passed  through  the  analyzer  exit  slit, 
they  entered  the  17  stage  secondary  electron  multiplier 
(SEM) .  In  a  SEM  a  positive  ion  is  further  accelerated  by 
the  potential  of  the  first  dynode.  On  striking  the  dynode, 
secondary  electrons  are  ejected  and  accelerated  to  the 
second  dynode  where  each  electron  ejects  more  electrons  etc. 
The  current  amplification  of  this  device  is  a  function  of 
the  voltage  across  its  dynodes.  The  approximate  absolute 
gain  is  given  in  Figure  2.7.  The  gain  was  measured  by 
collecting  the  positive  ion  current  on  the  first  dynode  of 
the  SEM  by  suitably  changing  the  dynode  potentials,  i.e., 
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Figure  2.7  Absolute  Gain  of  Secondary  Electron  Multiplier 
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the  first  dynode  was  at  earth  potential  through  its 
connection  to  the  ammeter.  The  second  dynode  had  a  small 
positive  bias  to  prevent  the  positive  ions  from  passing 
from  the  first  dynode  to  the  second. 

Normally  there  is  a  uniform  potential  drop  between 
each  dynode,  the  last  dynode  being  at  earth  potential  by 
its  connection  to  the  ammeter.  The  ratio  of  the  electron 
current  from  the  17th  dynode  of  the  SEM  to  the  positive  ion 
current  collected  from  the  first  dynode  is  the  absolute 
gain  of  the  SEM.  The  currents  were  measured  by  a 
Victoreen  Femtometer  with  a  10^  ohm  input  resistor. 

2 . 7  Gas  Handling  Plant 

Introduction  of  gas  into  the  ion  source  of  the  mass 
spectrometer  was  done  by  the  gas  handling  plant  (GHP) . 

Figure  2.8  shows  the  glass  reservoirs,  valves,  mercury 
manometers,  and  leaks  which  made  up  the  GHP.  Gases  were 
admitted  through  a  12/30  S?  ground  glass  connection.  The 
valves  were  a  combination  of  teflon  gasketed  metal  valves 
and  vacuum  cup  glass  stopcocks.  Dow  Corning  grease  was 
used  on  the  glass  stopcocks.  The  system  was  pumped  through 
a  cold  trap  by  a  lOL/sec  mercury  diffusion  pump  backed  by 
a  Welch  Duo-Seal  forepump. 

The  handling  of  water  and  ammonia  was  done  differently. 
The  ammonia  was  admitted  to  the  GHP  from  a  tank  of  anhydrous 
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Figure  2.8  Gas  Handling  Plant  for  introduction  of 


gases  into  the  ion  source. 
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ammonia.  The  gas  was  stored  in  the  2L  and  11L  reservoirs 
and  in  a  1L  bulb  which  was  connected  to  the  GHP .  Air  was 
removed  by  condensing  the  ammonia  in  the  cold  finger  of 
the  1L  bulb  and  pumping  off  the  non- condensable  gases. 

The  ammonia  was  admitted  to  the  ion  source  of  the 
mass  spectrometer  through  the  150  Lusec  leak* *.  The  leak 
not  only  drops  the  reservoir  pressure  but  also  prevents  a 
discharge  from  the  ion  source  (+1850  v)  to  a  grounded 
part  of  the  GHP.  The  pressure  in  the  ion  source  was 
varied  by  changing  the  ammonia  reservoir  pressure  by 
condensation  of  ammonia  in  the  cold  finger  of  the  1L  bulb. 

In  handling  water  vapor,  a  gas  bulb  filled  with 
about  50  ml  of  degassed  doubly- distilled  water  (liquid)  was 
used  as  a  constant  pressure  source  of  water  vapor. 
Temperature  control  of  the  bulb  by  a  dewar  filled  with 
water  or  ice-salt  mixture  varied  the  water  vapor  pressure. 
The  vapor  was  admitted  to  the  ion  source  by  the  70-100y 
sintered  glass  leak  (a  high  conductance  leak) . 

2 . 8  Measurement  of  Ion  Source  Pressures 

In  earlier  work  a  McLeod  gauge  had  been  used  to 
measure  the  ion  source  pressure.  Since  McLeod  gauges 
are  suitable  only  for  non-condensable  gases,  a  dibutyl 
phthalate  manometer  was  used  in  the  present  work.  Because 
there  was  a  pressure  drop  along  the  gas  inlet  line  between 

* 

One  Lusec  is  one  liter-micron/sec. 
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the  GHP  leak  and  the  ion  source,  the  dibutyl  phthalate 
manometer  had  a  separate  line  to  the  ion  source. 

When  the  present  work  was  started,  the  dibutyl 
phthalate  was  degassed  by  refluxing  it  under  vacuum.  To 
obtain  accurate  pressure  measurements  a  cathe tome ter  was 
used  to  read  the  height  difference  to  ±0.1  mm  which 
corresponds  to  ±0.0077  torr  ion  source  pressure. 

Because  ammonia  slowly  outgassed  from  the  dibutyl  phthalate 
in  the  manometer  am  under  vacuum,  the  vacuum  arm  was 
continuously  pumped. 

2 . 9  Temperature  Control  of  the  Ion  Source 

Ion  source  temperatures  above  room  temperature  were 
achieved  by  two  heaters  located  in  the  heater  holes  shown 
in  Figure  2.6.  Each  heater  was  made  by  winding  a  0.010 
inch  diameter  Mo  wire  around  a  quartz  rod  and  placing  the 
rod  and  wire  coil  inside  a  5  mm  diameter  quartz  tube. 

About  30W  was  needed  to  heat  the  ion  source  to  400°C  and 
60W  was  needed  to  reach  600°C. 

For  ion  source  temperatures  less  than  room  temperatures, 
two  different  cooling  methods  were  used.  During  the  low 

-f-  *f* 

temoerature  measurements  on  the  H  (HO)  and  NH  (NH  ) 

£  2  n  4  3  n 

systems,  the  ion  source  cooled  itself  by  a  net  radiation 
loss  to  the  liquid  nitrogen  baffles.  One  of  these 
baffles  was  the  one  leading  to  the  diffusion  pump  and 
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the  other  was  a  specially  made  baffle  consisting  of  a 
liquid  nitrogen  reservoir  attached  to  a  cylindrical 
copper  surface.  Together  the  baffles  virtually 
surrounded  the  ion  source.  Temperature  control  was 
possible  by  passing  only  a  very  small  current  through 
the  heaters. 

In  the  K  (H^O)^  stU(3ies  the  cooling  method  was  changed. 
The  heaters  were  removed  and  the  heater  holes  in  the  ion 
source  were  used  as  cooling  channels  of  a  circulatory 
system.  Two  Kovar  seals  were  soldered  to  the  two  heater 
holes  on  the  same  face  of  the  ion  source.  Each  of  these 
seals  were  glassblown  to  glass  tubing  which  led  out  of 
the  vacuum  housing  through  "0"  ring  seals. 

To  complete  the  ion  source  circulatory  system 
copper  tubing  was  connected  to  the  remaining  open  ends  of 
the  heater  holes.  Finally  the  system  was  checked  for 
vacuum  tightness.  Either  water  or  acetone  was  circulated 
through  the  ion  source  by  a  Colora  Ultra-Thermostat  which 
regulated  the  bath  temperature  to  ±0.01°C.  To  achieve 
temperatures  below  0°C,  the  bath  was  filled  with  acetone 
and  cooled  with  a  2L  stainless  steel  beaker  which  was 
filled  with  dry  ice.  The  depth  of  beaker  immersion  in 
the  bath  varied  the  temperature  to  as  low  as  -77°C.  The 
temperature  in  the  ion  source,  however,  went  only  to 


-50°C. 
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2 • Ion  Source  Temperature  Measurement 

The  temperature  of  the  ion  source  was  measured  by  iron- 
constantan  thermocouples  made  from  number  30  A.W.G.  wire. 

The  small  diameter  wires,  0.010  inch,  were  used  to  prevent 
an  erroneous  reading  due  to  heat  conduction  along  the 
wire.  The  locations  of  hie  thermocouples  are  given  in 
Figures  2.4  and  2.6.  In  the  ion  source  used  with  the 
proton  beam,  the  thermocouples  were  placed  in  wells  as 
shown  in  Figure  2.6.  Since  the  distance  from  the  bottom 
of  each  well  to  the  ion  source  wall  whose  temperature 
was  desired  to  be  measured  was  about  0.5  mm,  the  arrangement 

■j* 

was  changed  before  the  K  (HO)  study  was  begun.  To 

2  n 

obtain  a  more  accurate  measurement  of  the  ion  source 
temperature  at  the  wall,  the  existing  thermocouple  well 
in  the  ion  source  bottom  was  drilled  through  to  the 
ion  source  cavity  and  capped  by  a  piece  of  0.005  inch 
thick  nichrome  sheet  which  was  spot  welded  to  the  wall. 

A  second  thermocouple  well  as  shown  by  Figure  2.4  was 
made  in  the  same  manner. 

The  thermal  emf  was  measured  relative  to  a  0°C 
reference  thermocouple  by  a  Honeywell  potentiometer.  The  » 

temperature  of  the  thermocouple  could  be  read  to  within 


1°C. 
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3-  ION  SAMPLING  IN  THE  HIGH  PRESSURE  MASS  SPECTROMETER 

In  this  research  project  the  thermodynamic  functions,  AG°  , 
Ah0 ,  and  AS°  for  gas  phase  solvation  reactions  were  measured.  This 
required  a  knowledge  of  the  relative  ionic  concentrations  in  the 
ion  source.  In  this  section  the  general  events  such  as  ionization, 
recombination,  and  diffusion  which  affect  the  ionic  concentrations 
are  described.  In  addition  the  possible  relative  changes  in  the 
ionic  concentrations  occurring  in  sampling  and  measurement  are 
discussed.  In  general  these  possible  changes  alter  the  ion 
ratios  and  therefore  affect  AG°  and  AS°  but  not  AH°. 

3. 1  General  Reactions  in  the  Ion  Source 

A.  Reactions  Initiated  by  the  Proton  Beam 

Thermal  energy  ions  for  the  study  of  NH  „+ (NH.J  and  H+(H^O) 

4  3  n  2  n 

equilibria  were  created  by  the  processes  given  in  Table  3.1.  The 
irradiated  gases  and  beam  energies  shown  in  the  table  are  not  the 
same  as  those  used  in  the  present  experiments  since  information 
pertinent  to  this  work  is  not  available  in  the  literature. 

However,  the  literature  data  are  sufficiently  representative  to 
provide  a  qualitative  picture  of  the  processes  to  be  expected 
in  the  present  system  in  which  water  or  ammonia  vapor  was 
irradiated  by  100  kev  protons.  Since  each  100  kev  particle 
undergoes  more  than  one  collision  at  the  high  pressures  used, 
reactions  3.3,  3.4,  and  3.5  occur  in  addition  to  reactions 


3. 1  and  3.2. 
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TABLE  3.1 


Reactions  of  a  High  Energy  Beam 


Reaction* 


+ 

3.1 

H 

+ 

N2 

— - -> 

[N  ] 

+ 

+ 

3.  2 

H 
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V 

- y 

[N2  1 

4* 

3.3 

H 

+ 

°2 

- y 

°2  1 

-y 

3.4 

H 

+ 

°2 

- ». 

l°2  J 

-y 

3.5 

H_ 

+ 

Ar 

- y 

Ar  + 

Energy 


-  kev 

•y 

+ 

+  H  +  e 

100 

-y 
+  H 

100 

-y 

+  H~ 

30 

~y 

+ 

+  H 

30 

:  +  e 

40 

Total 

Cross 0Section  Refer- 

A  2  ence 

6  1 

3  1 

0.3  4g 

3  41 

20  42 


The  bracketed  products  represent  the  sum  of  all 

+  +  +  + 

the  low  energy  ion  products  l . e .  ,  N  ,  ,  O  , 
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The  thermal  energy  ions  created  by  the  high  energy  beam 
diffuse  to  the  walls  of  the  ion  source  where  they  are  discharged. 

En  route  to  the  wall  the  ions  react  to  form  solvated  hydrogen 
ions.  Some  of  the  ions  undergo  recombination  with  the  electrons 
or  negative  ions  created  by  the  processes  represented  by  reactions 
3. 1,  3.4,  and  3.5. 

To  calculate  the  approximate  steady  state  concentration 

of  the  positive  ions,  it  is  necessary  to  know  the  composition 

of  the  beam  and  the  cross  sections  for  reaction.  Since  these 

facts  are  not  well  known,  the  approximate  steady  state  concentration 

of  the  thermal  energy  ions  cannot  be  calculated  with  any  accuracy. 

12 

One  can  estimate,  however,  that  about  10  positive  ions/sec 

are  created  at  1  torr  pressure,  298°K  in  the  plane  of  the 

high  energy  beam  directly  above  the  ion  exit  slit  whose 

length  is  4  mm.  The  ion  production  calculation  assumes  a 

°2 

cross  section  of  3A  and  uses  a  beam  intensity  corresponding 
to  the  proton  beam  current,  0.05  yA,  measured  in  the  Faraday 
cup  of  the  ion  source  at  background  pressure.  Under  these 
conditions,  each  high  energy  particle  undergoes  about  10 
collisions  per  centimeter  path  length. 

Since  the  steady  state  concentration  of  positively 
and  negatively  charged  particles  may  be  high  in  the 
ionization  region,  recombination  may  be  a  major  loss  of  ions 
in  this  region.  Hex’/ever,  since  there  is  a  long  path  length, 

0.5  cm  to  the  ion  exit  slit,  recombination  in  the  narrow 
ionization  region  should  not  affect  the  relative  ionic 


' 
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cluster  concentrations  at  the  exit  slit. 


B.  Production  of  the  Solvated  Potassium  Ion 
To  produce  potassium  ions  for  the  study  of  hydrated 
potassium  ions,  a  heated  filament .was  utilized  rather  them 
the  proton  beam.  Consequently  the  initial  reactions  are  different. 
The  potassium  ions  after  having  been  ejected  from  the  heated 
filament  were  attracted  to  the  ion  source  by  an  electric  field 
of  approximately  lOOv/cm.  The  high  field  was  required  to 
attract  a  sufficient  flow  of  ions  into  the  ion  source  to  allow 
ion  current  measurement  after  mass  analysis.  Upon  entering 
the  field  free  ion  source  the  potassium  ions  had  excess 
translational  energy  i.e.  energy  greater  than  thermal  energy. 

Because  the  filament  chamber  operates  at  the  same 
pressure  as  the  ion  source,  the  translational  energy  of  the 
potassium  ion  can  be  calculated  by  Wannier's  equation  ( 43  ) , 


3.1 


E. 

ion 


1 

2 


(m.  +  m  ,  ) 

ion  neut. 


2  3 

V  +  /2  kT, 


for  the  translational  energy  of  an  ion  at  low  X/P  where  X 
is  the  electric  field  strength,  P  is  the  pressure,  and  V  is 

the  drift  velocity  of  the  ion. 

The  ion  drift  velocity  can  be  calculated  from  the 
mobility,  the  drift  velocity  per  atm  -  v/cm.  Munson,  Tyndall, 

■f 

and  Hoselitz  (  44)  have  measured  the  mobility  of  K  in  the 
gases  with  up  to  2. 8%  water  vapor  added  at  X/P  up  to 
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30v/ cm  -  torr.  For  values  of  X/P  from  0  to  10  v/cm  -  torr, 
the  mobility  was  constant.  Above  lOv/cm  -  torr,  the  mobility 
increased  with  increasing  X/P.  The  effect  of  increasing 
the  percentage  water  vapor  was  the  extension  of  the  constant 
region  of  the  mobility.  Assuming  the  constant  region  extends 

to  100  v/cm  -  torr  at  100%  water  vapor  gives  a  value  of  the 

.  .  2  -1  -1  -1 

mobility  of  2.2  cm  sec  v  atm  by  extrapolating  the 

data  for  the  mobility  of  K+  in  moist  Ar.  Under  the 

experimental  conditions  of  100  v/cm  -  torr  the  drift 

5 

velocity  is  then  1.7  x  10  cm/sec.  On  substitution  of 
the  drift  velocity  value  into  equation  3.1,  the  potassium 
ion  energy  is  0.74  ev.  This  means  about  six  collisions  are 
required  in  the  ion  source  to  remove  the  excess  translational 
energy  assuming  each  collision  removes  half  of  the  excess 
energy. 

The  maximum  depth  to  which  an  ion  penetrates  into  the 
ion  source  can  be  calculated  from  the  mean  free  path.  The 
formula  for  the  mean  free  path  A  is 


3.2 


A  = 


(k  °A  +  *  V 


where  N  is  the  neutral  gas  density,  c  is  the  collision 
B  “ 

diameter  of  the  ion  and  a  the  collision  diameter  of  the  gas 

B 

o 

molecule  (  26  ) .  Taking  0  as  2.66A  for  the  ionic  diameter  of 

A 

.  ° 

K  (  )  and  a  as  2.72A  for  the  diameter  of  the  water 

molecule  (27  )  gives  A  equals  0.01  cm  at  1  torr,  298°K.  At 
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1  torr  the  penetration  of  the  unsolvated  ion  into  the  ion  source 
is  less  than  1  mm.  At  0.1  torr  the  penetration  increases  to  8  mm 
as  calculated  by  the  above  considerations.  These  distances  can 
be  compared  to  the  1.6  cm  distance  from  the  potassium  ion 
entrance  orifice  which  is  in  the  ion  source  top  to  the  ion 
exit  slit.  Thus,  even  at  the  lowest  pressures  used,  the  path 
length  at  which  the  ion  is  at  thermal  energy  is  3  mm  longer  than 
the  path  length  in  the  study  of  the  and  NH^+(NH^)n 

systems . 

3 . 2  Transport  of  the  Solvated  Ions  Out  of  the  Ion  Source 

After  the  potassium  ions  lose  their  excess  energy  and, 
in  the  case  of  the  ions  created  by  the  proton  beam,  after  the 
ions  are  created,  the  ions  react  with  neutral  molecules  as 
they  diffuse  toward  the  wall.  During  this  time  the  ions 
should  form  equilibrium  ionic  clusters.  The  average  distance 
traversed  by  an  ion  along  an  axis  X  is  given  by 

~~2 

3.3  AX  2Dt 

2 

where  AX  is  the  average  mean  square  displacement,  D  is  the 

diffusion  coefficient,  and  t  is  the  time  of  diffusion  (24  ) . 

2  2 

Setting  AX  =  X  as  an  approximation  gives  upon  rearranging 

2 

3.4  t  =  X_ 

2D 

Substitution  of  0.5  cm,  the  distance  between  the  proton  beam 
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plane  and  the  ion  exit  slit,  for  X  and  taking  50  cm ’/sec  as 
a  representative  diffusion  coefficient  at  1  torr,  298°K  (25  ) 
leads  to  a  diffusion  time  to  the  ion  exit  slit  of  2.5  millisec. 
This  time  is  also  an  estimate  of  the  ion  source  residence 
time  of  an  ion.  This  estimate  is  approximate  since  space 
charge  effects  due  to  presence  of  the  charged  particles 
may  alter  the  diffusion  coefficient. 

The  calculation  applies  to  a  static  gas.  In  the  ion 
source  the  gas  flows  out  the  ion  exit  slit.  The  flow  rate 
can  be  approximately  calculated  by  modifying  an  equation  for 
the  flow  out  of  a  circular  orifice  (  46  ) .  The  flow  rate 
is  derived  from 

3.5  F  =  -  A  dr 

dt 

where  F  is  the  conductance  of  the  slit,  A  is  the  area  of  the 
hemi cylindrical  surface  illustrated  in  Figure  3.1  and  dr/dt 
is  the  flow  velocity  at  a  distance  r  from  the  center  line  of 
the  ion  exit  slit.  On  substitution  for  the  area  and 
rearrangement,  one  obtains  equation  3.6 

3.6  dr  =  -  F 

dt  irlr 

where  1  is  the  length  of  the  slit.  From  equation  2.6  the 

3 

conductance  of  the  ion  exit  slit  is  2  cm  /sec  for  a  molecular 
mass  of  18,  298°K  and  slit  dimensions  of  3y  x  0.4  cm.  Use 
of  equation  3.6  indicates  that  although  the  gas  is  flowing, 
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Figure  3.1 


Assumed  Gas  Flow  in  the  Ion  Source. 
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the  rate  does  not  become  significant  (lOOcm/sec)  until  about 
1  mm  from  the  exit  slit.  The  ion  residence  time  is,  therefore, 
not  greatly  shortened  due  to  the  gas  flow. 

The  major  transport  of  ions  then  occurs  by  diffusion. 
Since  the  diffusion  coefficient  is  inversely  proportional 
to  the  molecular  mass,  ions  of  lower  mass  should  have  a 
lower  ion  residence  time.  This  increases  the  observed  per 
cent  total  ionizationof  the  low  mass  ions.  However,  if  the 
equilibrium  ionic  clusters  rapidly  gain  and  lose  ligands 
(i.e.  fast  forward  and  reverse  reaction  rates),  all  ions 
will  have  the  same  average  molecular  mass.  The  diffusion 
process  thus  need  not  give  a  discrimination  in  the  ion 
sampling  process. 

One  process  that  may  give  a  discrimination  is  the 
passage  of  ions  out  of  the  exit  slit.  By  equation  2.6  low 
mass  ions  are  favored  to  exit  the  ion  source.  However,  the 
high  gas  flow  rate  at  the  ion  exit  slit  may  tend  to  lessen 
this  discrimination.  By  equation  2.6  the  flow  rate  at  the 

4 

slit  is  1.5-x  10  cm/sec  for  molecular  mass  18,  298°K. 

3. 3  Collisions  Outside  the  Ion  Source 

After  the  ions  and  gas  molecules  have  emerged  from  the 
ion  source,  any  ion-molecule  collisions  will  lead  to  non¬ 
equilibrium.  This  problem  was  encountered  by  Collins  (22) 
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in  a  preliminary  study  of  the  equilibrium  reaction 

3.6  H+(H  0)  .  +  H  O  - - *  H**~  (H  O)  . 

2  n-1  .  2  -  2  n 

Collins  found  that  the  ion  cluster  size  began  to  decrease 
after  a  pressure  of  about  0.3  torr  due  to  ionic  cluster- 
water  molecule  collisions  outside  the  ion  source.  Dissoc¬ 
iation  occurred  if  a  collision  occurred  since  the  ionic 
cluster  was  electrically  accelerated  upon  leaving  the  ion  source. 
Since  a  larger  hydrate  has  weaker  ion-ligand  bond  strengths 
and  may  have  a  larger  collision  cross  section,  the  larger 
hydrates  were  preferentially  dissociated.  The  extraneous 
collisions  were  eliminated  by  increasing  the  pumping  speed 
at  the  ion  exit  slit.  The  pumping  speed  at  this  point  was 
limited  by  the  low  gas  conductance  caused  by  the  close 
spacing  of  the  ion  source  and  the  first  accelerating  electrode. 
To  improve  the  conductance,  the  first  accelerating  electrode 
was  moved  from  a  distance  of  1  mm  from  the  ion  source  to 
a  distance  of  10  mm.  Because  the  conductance  is  proportional 
to  the  distance  squared  for  this  pumping  geometry  ( 22  ) , 
the  conductance  was  100  times  greater.  Figure  3.2  compares 
Collins'  results  with  the  results  from  this  work. 

As  shown  in  Chapters  4,  5,  and  6,  the  calculated 
equilibrium  constants  did  not  change  with  pressure  up  to 
some  maximum  pressure  P  after  which  the  equilibrium  constant 
began  to  decrease.  Evidently  at  higher  pressures  dis¬ 
equilibrium  collisions  again  become  significant. 
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At  higher  temperatures  it  was  found  that  this  maximum 
pressure  P  increased.  This  may  be  partly  due  to  fewer 
molecules/sec  effusing  out  of  the  ion  source.  As  the 
temperature  is  increased,  the  conductance,  vol/sec 
increases  by  T  by  equation  2.6  but  the  molecules/sec 
effusing  out  of  the  ion  source  per  sec  decreases  by  T  ^  by 
the  ideal  gas  law.  In  addition,  the  cluster  sizes  are 
smaller  at  higher  temperatures  and  so  the  cross  sections 
are  expected  to  be  smaller. 

3.4  Ion  Discrimination  in  Acceleration,  Mass 

Analysis  and  Detection 

In  the  preceding  sections  some  of  the  possible 
pitfalls  in  the  ion  sampling  were  described.  In  addition 
the  mass  spectrometer  itself  can  alter  the  ion  ratios. 

During  the  ion  acceleration  either  high  or  low  mass  ions 
may  be  discriminated  against.  The  fact  that  two  mass 
spectrometers,  the  proton  beam  mass  spectrometer  (operated 
by  A.  Zolla  and  the  author)  and  the  alpha  particle  mass 
spectrometer  (20  ) (operated  by  M.  Arshadi  and  J.  Scarborough) , 
with  two  very  different  sets  of  acceleration  electrodes 
gave  similar  results  from  the  study  of  reaction  3.6  suggests 
that  the  discrimination  in  ion  acceleration  was  not  very 
important. 

Discrimination  may  also  occur  in  the  secondary 
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electron  multiplier.  The  number  of  electrons  ejected  from 
the  first  dynode  by  the  impacting  ion  can  vary  with  the 
ion  velocity.  To  measure  the  SEM  discrimination  the 
relative  gain  of  the  SEM  was  measured  for  two  series  of 
solvated  ions.  The  procedure  is  given  in  Section  2.6. 

From  these  measurements  the  relative  gains  for 
H+(H20),  H+(H20)2,  H+(H20)3,  and  H+(H20)4were  (1.0), 

1.1,  0.8,  and  0.8  respectively.  For  NH^  (NH  )  , 

NH+(NH  )  ,  and  NH  +  (NH  )  the  relative  gains  were  (1.00), 

*4" 

0.95  and  0.95  respectively.  In  the  case  of  the  K  (H  O) 

^  n 

series  the  ion  currents  were  too  small  to  measure  directly 
However  it  appears  that  the  SEM  introduces  a  maximum 


discrimination  of  20%  in  the  measurement  of  ionic  clusters 
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4  *  RESULTS  and  discussion  of  the  hydration  of  k+  and 

COMPARISON  OF  RESULTS  TO  ELECTROSTATIC  CALCULATIONS 


4.1  Introduction 

Potassium  ions  which  are  ejected  by  a  heated  filament 
drift  into  the  ion  source  under  the  influence  of  an  electric 
field.  In  the  field  free  ion  source  the  potassium  ions  lose 
their  excess  translational  energy  and  react  with  water 
molecules  by  the  following  reactions: 


4.1  (0,1) 

4.2  (1,2) 

4.3  (n-l,n) 


+  M  + 

K  +  H20  - K  (H20) 

+  M  + 

K  (H20)  +  H20  - »  K  (H20)2 


K+(H20)n-l  +  V 


M  + 

- ►  K  (HO) 

« -  2  n 


Since  the  attachment  of  a  water  molecule  to  an  ion  is  exothermic, 
third  bodies,  which  are  indicated  above  as  M,  are  required  for  the 
removal  of  excess  energy. 

The  equilibrium  constant  for  the  general  reaction  is  given 


by  equation  4.1  where  IK+(H  0)  Wy)^ 


4.1 


Kn-l,n  IK+(H20)n 


IK+(H2°)n.1PH20 


are  the  mass  analyzed  ion  currents.  The  ion  currents  are  assumed 

proportional  to  the  corresponding  equilibrium  concentrations  of 

K+ fH  0)  and  K+(H  0)  in  the  ion  source.  Determination  of 
2  n  2  n-1 

K  at  different  temperatures  is  used  to  obtain  AH°n_1  and 

n-l,n  ' 

AS*  .  . 

n-l,n 
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4 . 2  Effect  of  the  Heated  Filament  and 

.  .  + 

Initial  Translational  Energy  of  K 

To  prove  that  the  relative  ionic  concentrations  were  not 
affected  by  the  presence  of  the  hot  filament  whose  temperature 
ranged  from  about  200°  to  600°  above  the  ion  source  temperature, 
three  different  filament  arrangements  were  used.  These  are 
illustrated  in  Figures  2.4  and  2.5. 

In  Source  I  the  filament  was  placed  in  a  separate 
compartment  attached  to  the  ion  source  top.  The  location  of 
Source  I  is  such  that  the  emitted  ions  and  infrared  radiation 
were  incident  on  the  ion  exit  slit.  To  show  that  the  results 
did  not  depend  on  non-equilibrium  heating  of  the  ion  source 
by  the  filament,  experiments  were  done  with  two  different 
diameter  orifices  between  the  filament  compartment  and  the 
ion  source.  The  larger  orifice  which  was  covered  by  an  open 
grid  had  a  diameter  of  5  mm  while  the  smaller  orifice  was 
1.5  mm  in  diameter.  As  can  be  seen  from  Figure  4.9  the 
agreement  using  the  two  arrangements  was  satisfactory  at 
higher  temperatures.  At  the  lower  temperatures  at  which 
K2  3  was  measured,  the  two  sets  of  points  begin  to  disagree. 
Apparently  the  ionic  concentrations  and  water  vapor  do  not 
achieve  thermal  equilibrium  with  the  temperature  of  the 
ion  source  wall  in  the  5  mm  orifice  arrangement.  At  333°K, 
for  example,  the  5  mm  orifice  point  would  be  on  the  ^ 
line  if  the  gas  were  8°  hotter  than  the  ion  source  wall 
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where  the  temperature  was  measured.  The  agreement  at  higher 
temperatures  is  not  surprising  since  under  these  conditions  the 
filament  temperature  is  closer  to  the  temperature  of  the  ion 
source  walls. 

Because  the  orifice  diameter  affected  the  results,  a 
second  source  of  K+,  Source  II,  was  installed  to  further 
check  the  results.  In  this  source,  the  emitted  infrared 
radiation  passed  over  the  ion  exit  slit  in  a  plane  5  mm 
above  the  slit.  Ions  entered  the  ion  source  in  the  same 
plane  but  could  diffuse  out  of  the  plane  after  having 
entered  the  field  free  ion  source.  These  ions  then  had  at 
least  5  mm  to  equilibrate  before  passing  through  the  ion 
exit  slit.  Figure  4.9  shows  that  the  results  obtained  with 
this  source  are  in  agreement  with  the  results  obtained  with 
the  small  orifice  in  Source  I. 

Another  way  in  which  the  effect  of  heat  transfer  to 
the  ion  source  can  be  examined  is  to  vary  the  filament 
wattage.  Usually  the  wattage  was  set  at  a  constant  value 
throughout  an  experiment.  The  minimum  wattage  was  used  to 
obtain  a  sufficient  ion  current  for  accurate  measurement. 

On  the  occasions  in  which  the  filament  wattage  was  changed 
during  an  experiment,  no  change  in  the  equilibrium  constant 
was  noted.  For  example,  at  58°C,  0.37  torr  the  wattage 
was  reduced  from  22W  to  14W.  ^  was  measured  at  both 
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wattages.  The  change  in  K  was  about  3%  which  is  approximately 
the  experimental  error  in  measuring  the  ion  currents. 

When  the  filament  wattage  was  changed  from  about  25W  to 
15W  in  a  study  of  the  solvation  of  K+  by  NH^,  the  equilibrium 
constants  decreased  by  about  30%.  For  this  reason  Source  III 
was  installed.  This  source  was  similar  to  Source  I  but  the 
filament  was  placed  much  farther  away  from  the  ion  source 
(see  Fig.  2.5).  The  use  of  Source  III  gave  results  in 
agreement  with  Sources  I  and  II  for  the  hydration  of  K+ 

i 

and  eliminated  the  change  of  equilibrium  constants  with 
filament  wattage  in  the  K+(NH^)^  sYstem*  The  agreement  of 
the  results  from  Sources  I,  II,  and  III,  therefore,  is 
evidence  that  the  heated  filament  does  not  introduce  non¬ 
equilibrium  thermal  conditions  in  the  ion  source. 

A  second  possibility  of  non- equilibrium  arises  from 
the  excess  translational  energy  with  which  the  potassium  ions 
enter  the  ion  source.  If  this  energy  is  incompletely  removed, 
the  measured  equilibrium  constants  will  be  too  low.  One  proof 
that  this  does  not  occur  is  the  agreement  of  the  results 
from  Sources  I  and  II.  In  Source  I  high  translational  energy 
ions  have  a  greater  chance  of  reaching  the  ion  exit  slit  and 
should  thus  be  sampled  with  greater  efficiency  than  low 
translational  energy  ions.  The  opposite  is  true  with  Source  II 
because  with  this  arrangement  the  ions  are  initially  directed 
over  the  slit  rather  than  at  it. 
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Figure  4.1  shows  a  plot  of  the  ion  ratio  K  (H^0)^/K  (H20) 

versus  the  electric  field  strength  used  in  Source  I  to  attract 

the  potassium  ions  into  the  ion  source.  The  invariance  of  '  . 

the  ion  ratio  plot  and  the  theoretical  calculations  for  the 

ion  translational  energy  in  Section  3.1B  prove  that  the 

equilibrium  constants  were  not  affected  by  ions  with  excess 

translational  energy. 

4. 3  Results  and  General  Discussion 
A.  Results 

Some  representative  data  are  shown  in  Table  4.1. 

The  mass  spectrum  is  tabulated  for  all  pressures  measured 
at  three  temperatures.  Figures  4.2  -  4.8  show  the  plots  of 

K  versus  water  pressure  obtained  at  different  temperatures. 

n-1, n . 

Van't  Hoff  type  plots  of  the  equilibrium  constants  are 
shown  in  Figure  4.9.  Table  4.2  lists  values  for  AH°29g,  Ag°298' 
and  AS°2g8  for  reaction  4.3.  The  values  for  AH0  and  AS°  were 
determined  from  equation  4.2  by  rreans  of  a  least  squares  treatment 

AH°  _  AS0  . 

4.2  log  K  =  n- 1  ,n  +  n-1  ,n 

n  ±,n  2.3RT  2.3R 


usinq  values  of  K  _  measured  at  the  various  pressures  and 
^  n-l/n 

temperatures.  In  the  base  of  reaction  5/6  where  there  is 
insufficient  data  to  apply  the  least  squares  treatment/  the 


thermodynamic  values  were  calculated  from  the  van  t  Hoff 

plot  shown  in  Figure  4.9.  AG°2g8  was  calculated  from  equation  4.3. 


, 
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CK+(H20)4D  /  ck+(h20)3d 
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Figure  4.1  Plot  of  K+ (H^O) ^/K+ (H^O) ^  with  Electric  Field 
Strength  Used  to  Attract  the  Potassium  Ions  into  the  Ion  Source. 


Electric  Field  Strength  v  /  cm 

Source  I  with  1.5  mm  orifice 
Pressure  0.37  torr 
Temperature  58° C 

Usual  Electric  Field  Employed:  lOOv/cm 
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Figure  4 . 3 
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Figure  4.4 
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P  (torn) 

Figure  4.5  Log  K  at  various  temperatures  and  pressures. 

f  -3 

The  data  (from  Source  I  with  a  5mm  orifice)  give  lower  free 
energy  values  than  the  data  in  Figure  4.4  due  to  non-equilibrium 


conditions . 
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Figure  4.6  Log  ^  at  various  temperatures  and  pressures 
£3  Source  I  with  1.5  mm  orifice;  0,  Source  II;  A 


/ 


/ 


0,  Source  II 


Source  III 


- 
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P  (torr) 


Figure  4.7  Log  ^  at  various  temperatures  and  pressures. 
L} ,  Source  I  with  1.5  mm  orifice;  Q ,  Source  II;  ,  Source  III 
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Figure  4.9  van't  Hoff  plots  of  equilibrium  constants  for  the 
hydration  of  K  :  0,  Source  I  with  5.0  mm  entrance  to  the  ion 

source;  O ,  Source  I  with  1.5  mm  orifice  entrance  to  the  ion 
source;  0,  Source  II;  and  A ,  Source  III. 
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TABLE  4 .  1 


Representative  Data  for  the  Hydration  of  K 


K+  Temp .  P . 

Source  C°  torr 

n=0 

Mass  Spectrum  of  K+(H^0) 

2  n 

in  %  Total  Ionization 

1  2  3  4  5 

6 

Total  Ion 

Current 

x  1016A 

P 

19  0.31 

42.5 

50.6  6.87 

0.23 

0.5 

0.51 

31.9 

54.4  13.2 

0.49 

0.2 

0.70 

24.7 

58.3  17.1 

0.07 

0.  84 

22.8 

56.4  20.4 

0.04 

II 

150  0.30 

66.4 

32.  7 

0.82 

0.1 

0.35 

65.2 

34.5 

1.28 

0.1 

0 . 66 

48.  7 

48.7 

2.66 

0.07 

0.88 

43.4 

52.8 

3.59 

0.05 

1.11 

36.7 

58.4 

4.87 

0.04 

1.  35 

35.4 

58.9 

5.59 

0.03 

ib 

264  0.42 

89.6 

10.2 

0. 10 

30 

1.28 

62.  3 

36.5 

1.20 

3 

2.56 

48.6 

48.6 

2.81 

0.8 

3.22 

39.1 

56.5 

4.33 

0.5 

3.64 

40.2 

55.0 

4.76 

0.3 

3.94 

35.2 

59.9 

4.94 

0.3 

5.15 

30.3 

69.7 

0.1 

3. 

Source  used  with  1. 

5  mm  orifice 

entrance  to 

the  ion 

source. 

b 

Source  used  with  5 

mm  open 

grid 

entrance  to 

the  ion 

source. 

. 
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TABLE  4.2 


Experimental  Thermodynamic  Values  for  the  Gas  Phase  Reactions: 


K* (HO)  + 
2  n-l 

H  O  — *  K+(H  0) 

2  ■< —  2  n 

(n-l,n) 

Reaction 

-AH°  , 
n«-l,n 

-AC°a 

n- 1 ,  n  ( 2-9  8  ° ) 

-AS  0  3 

n-l,n 

Ab 

Bb 

n-l,n 

(kcal/mole) 

(kcal/mole) 

e.u. 

i — 1 

o 

17.9 

11.4 

21.6 

3.  89 

7.57 

1/  2 

16.1 

8.9 

24.2 

3.50 

8. 14 

2,3 

13.2 

6.3 

23.0 

2.87 

7.88 

3,4 

11.8 

4.4 

24.7 

2.57 

8.26 

4,5 

10.7 

3.2 

25.2 

2.33 

8.  35 

5,6 

10.0 

2.3 

25.  7 

2.18 

8.46 

a  At  standard  state  of  water,  1  atm. 

b>  3 

Constants  for  the  equation  log  (K  )  =  (A/T)  x  10  -B. 

n  x  /  n 

Data  for  water  pressure  in  torr. 
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4.3  A.G°  =  AH0  -  TAS° 

B.  Entropy  of  the  K  (H  O)  Complex 
The  experimental  entropy  value  AS0  can  be  compared  to  a 

U  z  1 

theoretical  value  obtained  from  calculation  of  the  translational, 

rotational  and  vibrational  entropy  changes .  The  values  for  the 

other  AS°  .  could  also  be  calculated  but  a  large  number  of 
n-l,n  • 

approximations  would  have  to  be  made. 

The  translational  entropy  of  a  species  in  the  gas  phase 
is  given  by  the  Sackur-Tetrode  equation,  (  47)  4.4, 

4.4  Sc  =  3  R  In  M  +  5  R  In  T  -  2.311 

tr  2  2 

where  M  is  the  ordinary  molecular  weight  (g/mole) .  The  gas 
constant  R  must  be  substituted  in  cal/mole“ degree  to  obtain 
the  entropy  in  cal/mole-deg.  The  change  in  translational 
entropy  for  the  condensation  reactions  considered  here  is 
given  by  4.5. 

M  • + 

4.5  AS0  =  3  R  In  K  (H2°}  -  5  R  In  T  +  2.311 

tr  2  V  x  mh20  2 

On  substitution  for  reaction  0,1  at  298°K,  AS°^_  =  -  33.6  e.u. 

The  rotational  entropy  change  can  be  calculated  from 
the  rotational  entropy  of  H^O  and  K  (H^O)  .  The  appropriate 
equation  (  47  )  is  4.6, 

a  f.  c°  =  RTdflnQ  )  +  R  In  Q 

&  rot  rot  ^rot 


dT 
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where  Q  ,  the  rotational  partition  function  (48),  is  given  by 


4.7  Q_.,_  =  8tt  kT  (8tt2  kT  I  I  I  )h 


•rot 


x  y  z 


ah 

4.7  where  1,1,  and  I  are  the  three  moments  of  inertia  about 
x  y  z 

the  three  principal  axes  and  a  is  the  symmetry  number.  Differ¬ 
entiation  of  Qrot  with  respect  to  T  and  substitution  into  4.6 
gives  4.8. 


4.8 


S° 


rot 


=  3  R  +  R  In  Q 


rot 


For  the  rotational  entropy  of  water,  substitution  of  a  =  2 ,  and 

-40  2 

1^,  1^,  and  1^  =  1.91,  1.02,  and  2.94  x  10  g  cm  respectively 

leads  to  S°  =  10.1  e.u.  at  298°K. 
rot 

In  order  to  calculate  S°  ■  for  K+(H-0)  the  K+  -  0H„ 

rot  2  2 

bond  length  must  be  known.  As  shown  in  Section  4.4  a  reasonable 
distance  between  the  centers  of  the  potassium  ion  and  oxygen 

O 

atom  is  2.6A.  The  moments  of  inertia  are  then  I  =  1.91, 

-40  2 

I  =  147  and  I  =  149  x  10  g  cm  .  The  three  principal 

Y  Z 

-f* 

axes  for  the  HO  and  K  HO  are  shown  in  Figure  4.10. 

Z  ^  * 

Substitution  of  the  moment  of  inertia  into  4.8  leads  to 

S°  =  19.3  e.u.  with  a  =  2. 
rot 

The  vibrational  entropy  change  for  reaction  0,1 
can  only  be  approximately  calculated  because  the  fundamental 


vibration  frequencies  of  K  (HO)  have  to  be  estimated  to 


■ 
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(0.52,0.770) 


(0.52,-0.77,0) 


The  z  axis  is  mutually  perpendicular  to  the  x  and  y  axes  and 
passes  through  the  point  (0.04,  0,0). 


Y 


-o 

(-0.81,0,0) 


(2.88,0.77,0) 


(2.38,-0.77  0) 


Principal  Axes  of  K  (H^O) 


The  Z  axis  is  mutually  perpendicular  to  the  X  and  Y  axes  and 
passes  through  the  point  (0.02,  0,0) 


Figure  4.10  Principal  Axes  of  11^0  and  k'  (H20)  used  to  calculate 


the  moments  of  inertia. 
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calculate  S0^^.  To  simplify  the  calculations  the  assumption 
is  made  that  the  vibrational  entropy  of  H^O  is  unchanged  by 
bonding  to  K+.  The  coordinated  H^O  is  treated  as  a  rigid 
molecule  which  has  a  vibration  along  the  potassium- oxygen  axis. 

A  value  for  the  force  constant  of  the  vibration  can  be  estimated 
from  the  electrostatic  calculations  given  in  the  next  section. 

Thus  from  the  calculated  potential  energy  plot  shown  in 
Figure  4.19  the  force  constant  k  for  the  vibration  can  be 
calculated  to  be  5.8  x  10^  dyne/cm  by  fitting  equation  4.9  (49  ) 
to  the  bottom  portion  of  the  P.E.  Curve.  From  equation  4.10  (  49  ) 
where  v  is  the  fundamental  vibration  frequency  and  y  is  the 


4.9 


A  P.E.  =  1  k  (R  -  R  ) 

T  eq. 


4.10 


v  =  1 


2  7T 


k 

y 


12  -1 

reduced  mass,  v  is  calculated  to  be  8.47  x  10  sec  after 

substitution  for  k  and  y  which  was  approximated  as 

2  3 

(18)  (39)/  (18  +  39)  (6.02  x  10  )  g/molecule.  From  =  ^hv 

and  equation  4.11  (  47  )  ,  =  2. 8  e,u*  298°K. 


4.11 


R 


vib 


£o\-  R  In  (1-exp (-eo/RT) ) 


exp  (eoAt’  '  1  l\kT‘ 


In  addition  to  the  vibration  along  the  potassium-oxygen 
Qxlg i  -y.^02-0  gxe  two  torsional  vibrations  of  the  water  molecule 
in  the  complex  about  its  y  and  z  axes.  The  fundamental  vibration 


I  '  I 
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frequency  can  be  calculated  from  equation  4.12  (48)  where  n 
is  the  symmetry  number,  A  P.E.  is  the  potential  energy  barrier 
for  rotation,  and  I  is  the  moment  of  inertia.  The  values  of 

4. 12  v  =  n  [A  P.E.\  1 

2  it  \  21  J 

A  P.E.  can  be  calculated  from  the  coordinates  of  the  molecule 
in  its  minimum  energy  position  (Figure  4.10)  and  in  its 
maximum  energy  position  (rotation  by  180°  about  the  y  and  z 
axes  from  the  minimum  position)  by  use  of  equations  4.17, 

a 

4.18  and  4.20.  The  resulting  values  are  29  and  34  kcal/mole 
for  A  P.E.  and  A  P.E.  respectively.  This  leads  to 

y  z 

v  =  1.59  x  10^  sec  ^  and  v  =  1.03  x  10  ^  sec  \ 

y  z 

E  =^hv  =  0.76  kcal/mole  and  E  =^hv  =  0.485  kcal/mole, 
y  y  z  z 

and  therefore,  by  equation  4.11:  S0^  =  1.6  e.u.  and 

S°  =  2.4  e.u. 
z 

The  total  entropy  change  can  now  be  calculated  by 

summing  the  individual  contributions.  The  total  calculated 

value  is  -17.6  e.u.  compared  to  the  experimental  value  of 

-21.6  e.u.  The  difference  between  the  two  values  suggests 

that  the  experimental  value  may  be  somewhat  smaller  than 

the  true  value  for  AS0  . 

U  / 1 

C.  Distribution  of  K+(H„0)  Clusters 
- - - — - — - — - -  z—— n - - 

as  a  Function  of  Water  Pressure 
The  relative  concentrations  of  species  K+(H2°^n  can 


- 
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be  calculated  at  a  given  water  pressure  P  by  using  equation  4.13. 

n 

4.13  log  [K  (HO)  ]  =  n  log  P  +  Y  log  K. 

2  n  L  1-1,1 

i=l 

The  equilibrium  constants  K.  .  at  a  given  temperature  are 

l  -L,i 

conveniently  obtained  from  the  van't  -Hoff  straight  lines, 
equation  4.14. 

4.14  log  K.  .  .  =  AxlOOO  -  B 

i-l,i  - 

T 

Value's  for  A  and  B  are  given  in  Table  4.2.  The  use  of  equation 
4.13  necessitates  knowledge  of  one  of  the  cluster  concentrations. 
,  setting  the  concentration  of  Kh  equal  to  unity  the  other  concen¬ 
trations  can  be  calculated.  Setting  the  sum  of  all 
concentrations  equal  to  100,  one  obtains  the  relative 
concentrations  in  per  cent  of  the  total.  Figure  4.11 
illustrates  the  relative  concentrations  as  a  function  of 
pressure  at  T  =  300°K.  The  plots  show  a  characteristic 

behavior  which  is  also  observed  in  plots  arising  from  the 
+ 

H  (HO)  data.  At  low  pressures  where  the  low  hydrates  are 
2  n 

present,  the  difference  between  the  stability  (free  energy) 
of  the  hydrates  is  large  and  the  concentrations  of  one  or 
two  hydrated  species  are  dominant.  At  higher  pressures 
where  the  higher  hydrates  become  stable,  differences 
between  the  free  energies  of  the  hydrates  are  small 
and  a  number  of  hydrated  species  are  present  simultaneously. 


By 


' 


Total  Ionization  of  K*(HpO), 
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r 


Figure  4.11 


Plot  of  Relative  Concentrations  of  K+(H  0)  with 

2  n 


Log  P  at  300°K. 
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The  effect  is  more  pronounced  in  the  H+(H  0)  results. 

2  n 

D.  Comparison  with  Earlier  work 

The  equilibrium  constant  K  has  been  obtained  in 

0/1 

earlier  vzork  by  Chupka  (  50  )  .  The  possibility  of  determining 

this  equilibrium  constant  occurred  to  Chupka  while  he  was 

measuring  the  dissociation  energies  of  gaseous  alkali  halide 

complex  ions.  To  produce  the  ions  K  and  K  (H^O)  he  passed 

-3 

water  vapor  at  a  pressure  of  2x10  torr  through  an  oven 
which  was  coated  with  KCl  and  heated  to  570°C.  The  resulting 
ions  from  the  thermal  ionization  and  equilibration  were  mass 
analyzed.  The  value  obtained  in  this  manner  was  K  =  1.25 
at  570°C. 


Unfortunately  the  present  determinations  of  K 

U  / 1 

extend  only  to  385 °C  so  that  a  direct  comparison  of  the 

K  values  can  not  be  made.  However  a  fairly  accurate 

0,1 

extrapolation  can  be  obtained  by  use  of  equation  4.14. 

The  extrapolated  value  for  K  obtained  in  this  manner 

is  0.93  (both  K  values  relate  to  a  standard  state  of 

0/1 

one  atm).  The  agreement  between  Chupka' s  and  the  present 

K  is  fortuitously  close  considering  the  widely  different 

0,1 

experimental. conditions  .  Also  the  Chupka  constant  was 
obtained  at  a  single  pressure  which  was -not  directly 


measured/ but  estimated. 


86- 


E.  Proton  Affinity  of  KOH 

One  interesting  result  from  the  present  work  is  that 
AH°o  can  be  used  to  calculate  a  value  for  the  proton  affinity 
of  KOH.  The  heat  of  formation  of  K  (H^O)  in  the  gas  phase  is 
given  by  equation  4.15. 

4.15  AH  (K+H  0)  =  AH  (K+)  +  AH  (HO)  +  AH 

^  3-  -L  L//-L 

The  values  for  the  gas  phase  heats  of  formation  of  the 
potassium  ion  and  water  are  AHf  (K+)  =  122  and  AH^(H  0)  =  -57 
kcal/mole  (  51)  respectively.  Substituting  into  4.15  one 
obtains  AHfK+(H20)  =  +47  kcal/mole.  From  the  value  of  AHf 
of  KOH(g),  -58/ and  AH^  of  H+/  367  kcal/mole ,  (  51),  the  proton 

affinity  of  KOH,  represented  by  reaction  4.4  and  equation  4.16 
can  be  calculated  to  be  262  kcal/mole. 

4.4  H+.  .  +  KOH.  - - >  K+(fI  0)  .  . 

(g)  (g)  2  (g) 

4.16  AH  =  AH  K+(H  0)  -  AH  H+  -  AH  KOH 

ir  •  a  •  r  /  it  r 

From  a  literature  search  it  appears  that  the  proton  affinity 
of  KOH  has  not  been  previously  determined.  In  comparison  the 
proton  affinities  of  H^O,  CH^OH,  an<^  are 

193  kcal/mole  respectively  (  51  ) •  In  this  series  the  increase 
in  proton  affinity  is  due  to  the  increasing  polarizability  of 
the  substituent  group  (H  <  CH3  <  C^)  .  This  makes  the  oxygen 
more  negative  in  the  protonated  species  and  increases  the  bond 


* 
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energy  between  the  proton  and  the  oxygen.  The  oxygen  in  KOH  is 
expected  to  have  a  high  negative  charge  because  the  compound 
is  ionic.  Consequently  the  proton  affinity  should  be  high. 

4 . 4  Comparison  of  Measured  Enthalpies  and  Calculated 

Electrostatic  Potential  Energy  Changes  for 

Hydration  of  K+ 

A.  Introduction 

Classical  electrostatic  energy  calculations  were 
made  for  the  potential  energy  changes  which  occur  when 

water  ligands  bond  to  K+.  A  P.E.  values  were  then 

n-l,n 

compared  to  experimental  values  of  AH°  .  at  298°K. 

n-l,n 

Actually  the  electrostatic  energy  corresponds  best 
to  AH  at  0°K  because  of  neglect  of  all  thermal  motion  in 
the  calculations.  The  comparison  at  298°K  assumed  ACp  =  0. 

In  spite  of  this  assumption,  which  involves  a  considerable 
approximation,  it  is  felt  that  meaningful  conclusions  can 
be  made  from  a  comparison  of  the  electrostatic  calculations 
with  the  experimental  enthalpy  values. 

B.  The  Electrostatic  Calculations 

The  electrostatic  potential  energy  of  an  ion 
surrounded  by  molecules  can  be  expressed  by  the  following 
terms:  ion-dipole  and  ion-polarizability  attraction  and 

dipole-dipole  and  ion-atom  repulsion.  A  point  charge  model 
which  was  used  in  these  calculations  is  shown  in  Figure  4.12. 
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O 

Q-  =  +  1.00 
1 


Figure  4.12 


K+(H  0)  model  used  in  the  theoretical  calculations. 
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The  charge  assignment  resulted  from  the  separation  of  the 

dipole  moment  of  water  onto  the  three  point  atoms.  The 

model  and  equations  are  similar  to  those  used  by  Doyne 

and  Caldwell  (  52  )  who  considered  the  H„0  (H„0)  system. 

3  2  n 

Equations  4.17  -  4. 21^ which  are  written  in  a  general 

form,  were  used  to  compute  the  P.E.  of  K^H^O)^  over  a 
•f* 

range  of  K  -  OH^  distances  by  means  of  an  IBM  360  Computer. 

The  program  used  is  given  in  the  Appendix. 


4.17  Ion-Dipole 

n  p 

EE=  l  l 

j=l  i=l 


(334)  (Q.  )  (Q.  .) 

ion  i/j 

R.  . 
i*3 


4.18  Ion-Polarizability 

n  p 

EPOL  =  -  l  l 

j=l  i=l 


(167)«W  ai,i 

R4.  . 

1/1 


4.19  Dipole-Dipole 

n  p  &-1  P 

EDIP  =1111 
L-2  i=l  j=l  k=l 


(334)  (gi,i)(ek,d 

R.  .  .  . 

1/3/ k#Jl 


4.20  Ion-Atom 

n  p 

EREP  =11 
j=l  i=l 


A.  . 
i/3 


12 

R  .  . 

i/3 


4.21  En  = 


EE  +  EPOL  +  EDIP  +  EREP 


■ 

• 
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Q.  is  the  ionic  charge  and  Q.  .is  the  charge  on  the 
ion  3  i/D 

i—-  atom  of  the  j—  molecule  with  p  atoms  per  molecule  and 

n  molecules  per  cluster.  R.  .  is  the  ion- atom  distance 

i  /  3 

th. 

while  R.  .  _  .  is  the  distance  between  the  i —  atom  of 

i/  3,k,£ 

th  th 

j —  molecule  and  the  k— -  atom  of  molecule  £.  a  is  the 
atomic  polarizability  and  A  is  a  semi-empirical  constant 
for  the  ion-atom  repulsion  whose  evaluation  will  be  described 
later.  The  numerical  factors  in  equations  4.17  -  4.21  are 
adjusted  to  give  the  energy  in  kcal/mole  for  Q  in  atomic 


°  °  3 

units ,  R  in  A  and  a  in  A 


In  the  case  of  4.20  the 


numerical  factor  is  part  of  A  which  has  the  units  A 

-1 


12 


(kcal/mole) 


The  a  values  used  were  0.60  for  oxygen  and 


0.44A'3  for  hydrogen  (  53  )  • 

The  computed  values  for  the  ion-dipole  and  ion- 

+ 

polarizability  potential  energies  of  K  (H^O)  over  a  range  of 

R  (the  distance  which  corresponds  to  the  distance  between 

3,1 

the  centers  of  the  potassium  ion  and  oxygen  atom)  from  2.2 
to  4.1  A  are  tabulated  in  Table  4.3.  The  numerical  values 
show  that  in  this  distance  interval  the  ion-dipole  energy 
is  much  larger  than  the  ion-polarizability  energy.  Furthermore 
the  predominance  of  the  ion— dipole  energy  over  the  ion 
polarizability  energy  increases  rapidly  over  this  range. 


For  the  higher  hydrates  many  structures  are  possible. 
The  assumed  structures  are  shown  in  Figure  4.13.  In  these 
structures  all  the  K  -  OH^  are  equal.  (This 


. 
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TABLE  4 . 3 

Potential  Energies  (kcal/mole)  for  K+  -  H^O  as  a  Function  of  R0 

' r  “  "  Z  '  j  i  1 


R3,l 

EE 

EPOL 

EE  & 

A° 

(ion-dipole) 

(ion-polarizability) 

EPOL 

2.20 

-24.04 

-6.37 

-30.41 

2.30 

-22.14 

-5.42 

-27.56 

2.40 

-20.45 

-4.64 

-25.09 

2.50 

-18.95 

-3.99 

-22.95 

2.60 

-17.61 

-3.46 

-21.07 

2.70 

-16.40 

-3.01 

-19.42 

2.80 

-15.32 

-2.64 

-17.95 

2.90 

-14.33 

-2.32 

-16.65 

3.00 

-13.44 

-2.05 

-15.49 

3. 10 

-12.63 

-1.81 

-14.44 

3.20 

-11.89 

-1.61 

-13.50 

3.  30 

-11.21 

-1.44 

-12.65 

3.40 

-10.59 

-1.29 

-11.88 

3.50 

-10.02 

-1.16 

-11.17 

3.60 

-  9.49 

-1.04 

-10.53 

3.70 

-  9.00 

-0.94 

-  9.94 

3.  80 

-  8.55 

-0.85 

-  9.40 

3.90 

-  8.13 

-0.77 

-  8.91 

4.00 

-  7.75 

-0.70 

-  8.45 

4. 10 

-  7.38 

-0.64 

-  8.03 
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Figure  4.13 


+ 

Assur.ed  K  (HO) 

2  n 


Structures 


A  Structures 


n-2,  linear 


B  Structures 


One  molecule  of  the  A  structure  has  been  rotated  by 
90°  to  obtain  the  B  structure. 


All  3  molecules  of  the  A  structure  have  been  rotated  by 
90°  to  obtain  the  B  structure. 


n=4,  tetrahedral 
"top"  view 


The  three  molecules  which  are  trigonal  to  each  other 
have  been  rotated  by  90°  from  the  A  structure  to  obtain 


structure  B . 
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Figure  4.13  Cont'd. 


A_£ttructures  B  Structures 


n=5/  bipyramidal 


All  molecules  but  the  top  one  have  been  rotated  by 
90°  to  obtain  the  B  structure  from  the  A  structure. 


n=6,  octagonal 
"top"  view 


All  molecules  but  the  "top* one  have  been  rotated  by 
90°  to  obtain  the  B  structure  from  the  A  structure. 


- 
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ineans  that  the  total  energy  change  calculated  from  equations 

4-V* 

4.17,  4.18  and  4.20  for  addition  of  the  n — water  is  the 

same  for  each  addition  of  water  at  a  constant  ion-ligand 

separation  distance) .  In  general  the  most  symmetrical 

ligand  arrangements  were  used  to  minimize  the  dipole-dipole 

repulsion's.  Two  different  geometries  were  used  for  each 

structure  to  examine  the  effect  of  rotation  of  the  ligands 

about  their  dipole  axes.  Structure  B  was  obtained  from  A 

by  rotation  of  some  of  the  water  molecules  by  90°.  The  two 

structure  groups  did  not  change  the  computed  dipole-dipole 

repulsion  within  the  precision  of  the  computation  i.e.  less 

than  ±  0.05  kcal/mole.  Figure  4.14  shows  the  total  dipole- 

dipole  repulsion  for  each  cluster  calculated  from  equation 

4.19  as  a  function  of  R  .  The  increase  in  dipole-dipole 

3  f  -L 

-f- 

repulsion  due  to  the  addition  of  H20  to  K  (H2°^n-1  ^or  n==^“6 

can  be  calculated  for  any  assumed  R3  1  distances  for  K+(H20)n_1 

and  K (H  O)  by  using  Figure  4.14.  Figure  4.15  illustrates 
2  n 

the  increase  in  dipole-dipole  repulsion  calculated  from  the 

O 

plots  in  Figure  4.14  at  R  =  2.6A  which  represents  an 

average  radius  expected  for  these  clusters.  It  was 

intuitively  expected  that  the  potential  energy  difference 

should  show  a  smooth/near  exponential  increase  with 

increasing  n— l,n.  The  actual  plot  shows  an  S  shape. 

•phis  can  be  explained  by  the  following.  Ihe  number  of 

.  th 

additional  ligand  repulsions  per  addition  of  the  n 


' 
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Figure  4.14  Plot  of  the  total  dipole-dipole  repulsion 

energy  due  to  the  reaction  K+  +  nH  0  - >  K+(H  0)  . 

2  2  n 

Total  Dipole-Dipole  Repulsion  (kcal/mole) 

-»•  ro  Co 


' 


aED!P„  1  ^  (kcal  /  mo!e) 
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n-1,  n 


Figure  4.15  Plot  of  the  change  in  dipole-dipole  energy 

for  the  reaction  K+(H  0)  +  HO  - >  K+(H^0)  where 

2  n- 1  2  2  n 


R_  =  2.6A 

3  ,n 
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ligand  is  n-1.  This  fact  suggests  the  smooth  exponential 
increase  in  AEDIP  with  increasing  n-l,n.  Table  4.4  shows 
the  distances  between  the  centers  of  the  oxygen  atoms  in 

O 

the  assumed  structures  at  R  =  2.6A.  The  values  in 

3,n 

parentheses  are  the  number  of  distinct  interactions  with 
this  separation  distance.  From  the  table  it  can  be  seen 
that  the  oxygen  atom  separation  decreases  with  increasing 
cluster  size  in  an  irregular  manner.  This  is  a  consequence 
of  the  different  point  group  structures  for  different  n. 

The  irregular  decrease  apparently  leads  to  the  shape  of 
the  plot  in  Figure  4.15. 

A  complete  calculation  of  the  P.E.  should  consider  in 
addition  to  the  electrostatic  forces  used  above,  the  ion- 
ligand  and  ligand- ligand  repulsion  energies  and  the  ion- 
ligand  and  ligand- ligand  London  dispersion  energies.  The 
latter  are  quite  small  (less  than  a  few  kcal/mole)  (52,54) 
and  can  be  neglected.  The  ligand- ligand  repulsion  energy  can 

be  expected  to  be  less  than  the  ion-ligand  repulsion  energy 

-12 

since  the  repulsion  varies  with  R  (54,55)  and  the 
ligand-ligand  distance  is  in  general  50%  greater  than  the 
ion-ligand  distance.  Since  the  ion-ligand  repulsion  energy 
is  a  minor  contribution  to  the  total  P.E.,  the  ligand-ligand 
repulsion  energy  can  be  neglected. 

To  calculate  the  ion-ligand  repulsion  energy  from 

equation  4.20  it  is  necessary  to  assign  values  to  Aq,  a 
constant  for  the  ion-oxygen  atom  repulsion  and  Ah,  a 
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TABLE  4.4 

Intermolecular  Oxygen  Atom  Separation  Distances  for  =  2.6A0 

-  —  g  j  n —  -  ■ 


n 

2 


Distance* 

A° 


5.20(1) 


3 

4 

5 

6 


5.07(3) 


4.26  (6) 


3.67(6) ,  5.07(3) ,  5.2(1) 


3.67(12) ,  5.20(3) 


The  numbers  in  parentheses  are  the  number  of  equal,, 
distinct  oxygen-oxygen  distances  in  the  cluster. 
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constant  for  the  ion-hydrogen  atom  repulsion.  Minimization 


of  the  total  P.E.  equations  for  K+(H  0)  with  respect  to  R 

2  1,3 

leads  to  equation  4.22  from  which  A^  can  be  calculated 


4.22  3  P.E.  =  -  l  (334)  ^ion1  +  4 


3  R.  . 
i/l 


i=l 


R2.  . 

1/1 


3 

I 

i=l 


(167)  (Q  .  )  (a.  .) 

_ ion  1,3 

R5.  .  ' 

1/1 


12 


A 


H 


R 


13 


1/1 


R 


_H_ 

13 


+ 


A. 


0 


2,1 


R 


13 


3,1 


=  0 


as  a  function  of  R.  assuming  each  value  of  R.  corresponds  to 


the  equilibrium  ion-ligand  distance.  The  terms  containing  A 


H 

may  be  expected  to  be  negligible  compared  with  the  A^  terms. 
Firstly  A^  is  expected  to  be  larger  than  A^  since  the  repulsion 
depends  on  the  number  of  electrons  associated  with  the  atom  (  55  ) 
Also,  the  R  and  R  distances  are  greater  than  R 

J-  f  -L  £  f  -L  g  _L 

Neglecting  the  A  containing  terms  in  equation  4.22  one  obtains 

H 

an  equation  of  the  general  form  4.23  which  gives  Aq  as  a 


4.23 


*0  - 


function  of  R  .  This  function  is  shown  in  Figure  4.16.  The 

d)  /  1 

-f*  t 

potential  energy  of  the  cluster  K  (HO)  is  given  by  equation 
4.21  which  depends  on  the  parameters  R^  1  and  Aq  as  indicated 
in  equation  4.24.  In  order  to  obtain  a  value  for  Aq  and  R^  ^ 


E1  =  F(V  R3,l> 


4.24 


‘ 


c 

(A  )(kcal  /  mole) 
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R3,1  ^ 


Figure  4.16  Plot  of  the  ion-oxygen  atom  repulsion  constant 
as  a  function  of  R3  j  ln  this  plot  represents  the  ion- 

atom  equilibrium  distance  as  a  variable. 
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at  equilibrium  we  equate  E  with  the  experimental  AH  , 

^  %  0  / 1 

which  leads  to  4.25.  Expressions  4.23  and  4.25  represent 


4.25 


iH0,l  -  F(V  R3,l> 


a  system  of  two  equations  with  two  unknowns  (Ao/R  ,  ) 

0  3 , 1  eq 

which  can  be  solved  graphically.  To  do  this  R  is 

3,1  eq 

varied.  The  values  of  R  with  the  corresponding  numbers 
for  A^  (obtained  by  4.23)  are  substituted  into  equation  4.24. 
This  procedure  leads  to  the  plot  E  versus  R  shown  in 

3  /  _L 

Figure  4.17.  The  value  of  R  corresponding  to  the 

j  /  1 

O 

potential  energy  E  =  AH  =  -  17.9  kcal  is  2.56A,  the 

i-  yJ  f  -L 

equilibrium  distance.  The  associated  value  for  A  from 
1  0 

+  5  0 12 

Figure  4.16  is  A^  =  3.3  x  10  A  /kcal/mole.  This  derived 

value  of  A  will  also  be  used  to  calculate  EREP  and  E  for 
O  n 

the  higher  clusters. 

E  is  shown  as  a  function  of  R  in  Figure  4.19.  Also 
JL  *3  f  X 

included  in  the  Figure  is  a  plot  for  E  .  Similar  plots  were 

b 

done  for  the  other  clusters.  The  energy  minima  in  these 

plots  corresponding  to  the  total  calculated  binding  energy 

of  the  clusters  relative  to  K+  and  H^O  molecules  at  infinite 

distance  is  E°  .  Thus  AH  ,  should  be  compared  with 

n  n-l,n 

E°  -E°  =AE°.  Unfortunately  the  E°  values  can  not  be 

n  n-1  n  n 

too  accurately  estimated  from  the  E^  curves  (see  Figure  4.19). 
This  leads  to  a  large  error  in  AE0^. 

To  more  accurately  calculate  AE°n/the  following 


- 


Total  RE.  (kcal  /  mole) 
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2.4  2.5  2.6  2.7 


R3,1  ^ 

Figure  4.17  Plot  of  the  total  P.E.  using  the  values  of 

A  derived  from  Figure  4.16.  R  in  this  plot  is  the 
O  4/1 


equilibrium  ion-atom  distance. 
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Figure  4.18  Plot  of  the  ion-atom  repulsion  energy  with  the 

5  °  12 

derived  value  for  A  ,  3.3  x  10  A  /  (kcal/mole) .  R  is  the 

U  j  /  J. 


ion-oxygen  atom  distance. 


- 
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Figure  4.19  Plot  of  the  potential  energy  evolved 

for  the  reactions  K +  H  O - >  K+(H  0)  and  K+  +  6  H  0  — ►  K+(H  0) 

2  2  2  2 

Potential  Energy,  E  (Real /mole) 


(V) 
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procedure  was  adopted.  The  equilibrium  ion- ligand  distance 

O 

for  all  clusters  was  taken  as  2.6A.  This  approximation  is 

satisfactory  since  the  equilibrium  distance  increases  only 
°  +  ‘  + 

about  0.04A  from  K  (HO)  to  K  (H„0)  .  The  AE°  calculated 

2  2  6  n 

O 

at  2.6A  are  shown  in  Figure  4.20  for  n=l  to  n=6. 

C.  Comparison  of  Electrostatic  Potential 

Energies  with  Previous  Calculations 

The  calculations  given  in  the  previous  section  are 

similar  to  those  of  Garrick  (28  )  who  used  a  point  dipole 

model  to  calculate  the  P.E.  of  a  number  of  hydrated  ions 
+ 

including  K  (H^O)^.  The  formulas  were  similar  except  for 
ion-molecule  repulsion.  Garrick  used  equation  4.26  which 

.  —  ( v—  1) 

4.26  EREP  =  n  X  r 

v-1 

was  adapted  from  Lennard- Jones  calculations  of  inert  gas 
ion- atom  repulsions.  v  was  a  force  index  taken  to  be  10; 

A  was  a  calculated  force  constant;  and  n  was  the  number 
of  attached  water  molecules. 

Garrick's  a  priori  calculations  for  the  hexahydrate 
gave  the  equilibrium  distance  between  the  point  ion  and 

O 

point  molecule  as  2.58A,  EREP  as  18  kcal/mole  and 
EE  +  EPOL  +  EDIP  as  102  kcal/mole/  and  the  total  P.E.  as 
84  kcal/mole.  In  comparison  the  values  from  this  work  are 
2.56A  for  the  ion-oxygen  atom  distance,  19  kcal/mole  for 


■ 

■ 
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EREP,  105  kcal/mole  for  (EE  +  EPOL  +  EDIP) ,  and  86  kcal/mole 
total  P.E..  Our  values  are  semi-empirical  since  the  repulsive 
term  was  fitted  so  that  the  total  P.E.  at  the  equilibrium 
internuclear  distance  equalled  the  experimental  value  of 
AH  .  The  agreement  between  the  calculations  suggests  that 
the  point  dipole  model,  which  can  be  termed  point  ion- 
point  molecule  model  is  about  as  good  an  approximation  as 
is  the  point  ion-point  atom  model  used  in  the  present  work. 

Besides  Garrick,  Laidler  (  54  )  has  performed  point 
dipole  P.E.  calculations  for  a  number  of  complex  ions 

-f-  *f* 

including  K  (HO) ^ .  Unfortunately  the  K  (F^O) ^  results 
were  presented  as  a  total  of  a  cycle  dealing  with  the  heat 
of  hydration  of  the  potassium  ion.  However  the  individual 

energy  contributions  to  the  cycle  are  given  for  the  ion 

+  ° 

Fb  (HO),.  Since  the  ionic  radius  of  Rb  is  only  0.16A 
2  4 

larger  than  the  radius  of  K+  (  45  )  ,  one  would  expect  the 

•f* 

electrostatic  energy  of  formation  of  Rb  (HO)  to  be 

slightly  smaller  than  the  energy  of  formation  of  K+(H20)4. 

The  total  P.E.  for  Rb(H  0)  +  by  Laidler's  calculations  is 

-61  kcal/mole  while  the  corresponding  energy  for  K+(H20)4 

by  the  present  calculations  is  -66  kcal/mole.  One  can 

also  compare  Laidler's  values  for  the  repulsion  energy 

of  Fb  (H  0)  +  (-7.6  kcal/mole),  and  Li (HO)  +  (-9.1  kcal/mole) 

2  4  2 

to  that  of  K+(H  0)  +  (-12  kcal/mole).  Since  both  Laidler 

2  4 

and  this  author  used  the  same  dependence  of  the  repulsion,  i.e. 


■ 


' 
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equation  4.20,  Laidler's  value  of  A^  (which  was  not  given  in  the 
paper)  is  either  somewhat  lower  or  R  somewhat  higher  than  that 
derived  in  this  work. 


D.  Discussion  of  Calculated  E°  and 


■n- 


Experimental  AH0 
— 1 - - — — - n- 1 ,  n 

The  purpose  of  the  theoretical  calculations  was  to 

provide  an  explanation  for  the  observed  changes  in  AH^  ^ 

Figure  4.20  shows  a  plot  of  AH°  and  AE°  in  function 

n-1, n  n 

of  n.  According' to  the  theoretical  model,  the  decrease  of 
-  AE0^  with  n  is  due  to  an  increase  in  dipole-dipole 
repulsions.  As  can  be  seen  from  Figure  4.20  the  electrostatic 
calculations  result  in  values  of  AE0^  which  qualitatively 
agree  with  the  experimental  AH°^_^  Thus  the  changes  in 

AH0  ,  should  be  mostly  due  to  increases  in  dipole- 

n-l,n 

dipole  repulsions. 

In  comparing  the  two  plots  shown  in  Figure  4.20  more 
closely,  it  is  apparent  that  -AH°n_1  ^  decreases  more  rapidly 
initially  than  -AE0^.  This  effect  may  be  partly  due  to  the 
naive  model  used  in  the  calculations.  For  example  it  has 
been  proposed  that  the  lone-pair  oxygen  electrons  can  screen 
the  cation  charge  (  57  ) .  The  coordination  of  each  succeeding 
water  molecule  will,  therefore,  involve  a  smaller  release 
of  energy.  The  inclusion  of  this  effect  in  the  electro¬ 
static  calculations  would  be  expected  to  make  the  plot  of 


' 


kcai  /  moie 
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n  - 1 ,  n 

Figure  4.20  Comparison  of  the  theoretical  potential  energy 

o 

changes  to  the  experimental  AH  for  the  reaction 

n  -L  f  n 


K+(H20)n-l  +  H2°  ~ *  K"(H20,n 


-109- 


AE°  vs.  n-l,n  more  similar  in  shape  to  AH0  ,  vs.  n-l.n. 
n  n-l,n 

At  n-l,n  values  greater  them  5,6  it  is  expected  that 

AE°^  will  decrease  much  more  rapidly  due  to  rapid  increases 

in  dipole-dipole  repulsion.  As  stated  earlier  the  overall 

shape  of  AEDIP  vs.  n-l,n  should  be  an  exponentially 

increasing  function.  Since,  by  analogy  to  the  h"*  (H„0) 

2  n 

results  up  to  n  =  8,  the  AH°  ,  values  are  expected  to  be 

n-l,n 

almost  constant  at  high  n-l,n,  the  present  model  would  be 
unsatisfactory  at  high  n-l,n.  However,  under  these  conditions 
the  model  should  be  modified  to  allow  for  the  formation  of  an 
outer  shell  i.e.  one  in  which  the  incoming  water  molecule 
hydrogen  bonds  to  an  existing  ligand.  The  value  of  AE0^ 
can  be  calculated  from  the  existing  equations  and  use  of  the 
model  shown  in  Figure  4.21.  The  model  was  chosen  somewhat 
arbitrarily.  The  outer  water  molecule  has  been  placed 
at  a  distance  and  a  position  corresponding  to  a  normal 
hydrogen  bond  (58).  From  equations  4.17,  4.18,  4.20  and 
4.21  the  total  energy  release  in  coordination  of  the  outer 
water  molecule  is  -10.3  kcal/mole  of  which  -4.8  kcal/mole 
is  due  to  the  hydrogen  bond  formation.  The  value  -4.8  kcal/ 
mole  for  the  hydrogen  bond  formation  compares  to  an 
experimental  value  of  -5.0  kcal/mole  (  59  )  •  The  formation 
of  an  outer  complex  would  not  be  energetically  expected 
to  occur  until  equation  4.27  holds. 

4  27  -5.5  kcal/mole  +  AEDIP  £  -17.6  kcal/mole  +  AEDIP. 

outer  inner 


♦ 


■  .  . 
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K*  “O 


Figure  4.21  Model  used  to  calculate  the  potential 
energy  evolved  for  the  reaction 


K  (HO)  .  +  HO 

2  inner  2 


K  (H  O)  .  (H  O)  4 

2  inner  2  outer 
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In  4.27  AeDIP  is  the  increase  in  dipole-dipole 

outer  x  c 

repulsion  in  formation  of  the  hydrogen  bonded  ligand  while 

AEDIP.  is  the  increase  in  dipole-dipole  repulsion  in 

inner  c  c  * 

coordination  of  the  water  molecule  to  the  ion  at  the  usual 

O 

distance,  R_  =  2.6A.  The  value  -17.6  is  the  energy 
3  ,n 

release  exclusive  of  AEDIP  in  formation  of  the  ion-water 

O 

bond  at  2.6A  while  -5.5  is  the  energy  release  exclusive  of 

AEDIP  in  formation  of  the  outer  complex.  From  the  model  and 

the  usual  equations,  the  formation  of  the  outer  complex  is 

1.2  kcal/mole  energetically  unfavourable  relative  to  the 

til 

inner  coordination  for  the  addition  of  the  6  water 

til 

molecule.  This  suggests  that  the  addition  of  the  7 
water  may  result  in  an  outer  complex.  The  plot  of  AE0^  vs. 
n-l,n  would  then  be  fairly  constant  at  n-l,n  greater  than 
5,6  in  agreement  with  the  expected  trend  in  AH°n_^ 


- 
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5.  HYDRATION  OF  THE  PROTON 


SUMMARY 


+ 

The  hydrates  H  (HO)  were  observed  at  water  vapor 

pressures  from  0.1  to  6  torr  and  temperatures  of  25°  to  600°C. 

The  values  of  K  ,  AH0  ,  ,  AG°  _  ,  and  AS0 

n-l,n  n-l,n  n-l,n  n-l,n 

were  determined  for  the  reactions  H+(PI_0)  ,  +  H^O - >  H+(H  0) 

2  n-1  2  ~  2  n 

where  n  varied  from  2  to  7. 

These  results  are  in  general  agreement  with  the  values 
measured  on  a  second  mass  spectrometer  (AMS)  by  M.  Arshadi 
and  J.  Scarborough  (  20  )  .  The  AMS  w7as  operated  at  somewhat 
higher  pressures  over  the  temperature  range  from  -17°  to 
125 °C.  Consequently  the  results  from  the  AMS  were  for  the 


reactions  with  n  equals  4  to  8. 


■f* 

5 . 1  Reactions  Leading  to  H  (H^O)^ 

The  main  primary  ions  formed  in  the  100  kev  proton 

,|.  -f- 

beam  are  H  0  and  OH  .  These  ions  react  rapidly  to  give  HO  . 

2  d 


+ 

+ 

1 — 1 

• 

in 

«2°  +  H2° 

— ”  h3° 

+ 

OH 

5.2 

0H+  +  H20 

— "H30+ 

+ 

0 

Using  the  thermal  rate  constant  6.0  x  10  cm  /molecule  sec  (  60  ) , 
one  can  calculate  that  99%  of  the  primary  ions  react  at  1  torr 
in  0.3  ysec.  to  yield  H30+.  The  further  reactions  of  H30+  are 


reactions  5.3  -  5.7. 


. 
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+  M  + 

5-3  H30  +  H20  ►  H  (H20)2  (1/2) 

+  M  4- 

5.4  H  (H20)2  +  H20  H  (H20)3  (2,3) 

+  M  + 

5.5  H  (H20)3  +  H20  - *  H  (H20)4  (3,4) 

5.6  H+(H20)4  +  H20  h+(h20)5  (4,5) 

5.7  H+(HO)  .  +  H  0  -M->H+(H  0)  (n-1  ,n) 

2  n- 1  2  < —  2  n 


Since  the  attachment  of  water  molecules  to  the  ion  is 
exothermic,  third  bodies  indicated  above  as  M  are  required 
for  removal  of  the  excess  energy.  A  somewhat  more  detailed 
mechanism  is  given  by  equations  5.8  -  5.11 

5.8  Association  H+(H^O)  .  +  Ho0  - >  H+(H^0)* 

2  n-1  2  * - -  2  n 

5.9  Deactivation  H+(H  0)  *  +  Ho0  - ->  H+ (IT  0)  +  H_0* 

2  n  2  ■< - -  2  n  2 

5.10  Dissociation  H+(Ho0)  +  Ho0  H+(Ho0)  .  +  2H  0 

2  n  2  < -  2  n-1  2 

5.11  Dissociation  H+(H  0)  +  H.O  H+(Ho0)  .  +  (H_0)  _ 

2  n  2  < -  2  n-1  2  2 

H+(H20)  *  represents  the  species  still  containing  the  energy 
released  by  the  condensation  reaction.  Dissociation  may  occur 
by  the  reverse  reaction  of  5.8,  by  5.10  or  by  a  lcwer  energy 
path  (relative  to  5.10)  through  hydrogen  bond  formation  as 
represented  by  reaction  5.11. 

Provided  that  equilibrium  is  achieved,  the  relationship 
given  by  the  following  equation  should  be  observed. 
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5.1 


K  —  X  + 

n-l,n  H  (HO) 

2  n 


V(H20)n-lPH2° 


H  (HO)  H  (HO)  are  the  mass  analyzed  ion  currents 

2  n  2  n-1  J 


•f*  -f- 

of  H  (H  O)  and  H  (H  O)  respectively.  P  ^  is  the  water 
2  n  2  n-1  J  H20 


vapor  pressure  in  the  ion  source. 


5.2 


Results 


A.  Values  for  K 
- ■  -  •  ■ - n- 1 ,  n 

Some  representative  data  used  to  calculate  K  .  is 

n-l,n 

shown  in  Table  5.1.  The  mass  spectrum  is  tabulated  for  various 
pressures  and  temperatures.  Figures  5.1  -  5.7  show  all  the 
measured  values  for  the  equilibrium  constants  as  a  function  of 
pressure.  As  illustrated  by  the  plots,  the  equilibrium  constants 
do  not  change  significantly  over  the  experimental  pressure  range 


B.  Values  for  AH0  ,  AG°  , 

- — - - - — n-1  ,n  — n-1  ,n 

and  AS°  ,  at  298°K 
- n- 1 ,  n - - - - 

Values  for  AH°  ,  were  determined  from  van't  Hoff 

n-1  ,n 

plots,  log  K  ,  versus  1/T.  The  plots  are  shown  in  Figure 
^  J  n-l,n 

5.8.  The  results  obtained  by  M.  Arshadi  and  J.  Scarborough 
(  20  )  ,  on  a  different  mass  spectrometer  (AMS)  are  also  shown. 
According  to  the  van't  Hoff  relationship,  equation  5.2,  the 
slope  of  the  plotted  lines  equals  -AH0  /2.3R.  Since 
within  experimental  error  AH0  was  constant  with  temperature. 


' 
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TABLE  5.1 


No. 

Temp 

a 

Representative  1 

Data  Given  as  the 

Mass  i 

Spectra 

Total  Ion 

Current 

A  x  1014 

.  *  C° 
b 

in  Per 

p. 

torr  n=1 

Cent 

Mass 

2 

Total  Ionization 

.  -f" 

Spectrum  of  H  (H^ 

3  4  5 

0) 

n 

6 

7 

1 

55 

63 

0.19 

3.15 

81.3 

15.3 

0.29 

4 

0.30 

74.3 

25.0 

0.70 

0.015 

7 

0.53 

65.2 

33.2 

1.71 

0.029 

8 

0.59 

58.4 

39.3 

2.37 

0.036 

6 

0.76 

54.2 

42.6 

3.32 

0.041 

4 

1.50 

46.4 

48.6 

5.0 

0.081 

0.6 

55 

63 

2.10 

42.3 

51.4 

6.4 

0.2 

2 

93 

105 

0.18 

31.6 

67.5 

1.15 

7 

0.34 

16.8 

80.4 

2.77 

12 

0.47 

13.  3 

81.8 

4.92 

18 

92 

105 

0 . 65 

11.9 

83,0 

5.16 

6 

0.92 

7.67 

83.6 

8.66 

6 

1.45 

8.92 

81.0 

10.0 

2 

1.60 

9.65 

78.8 

11.6 

1 

93 

108 

2.30 

10.3 

75.5 

14.3 

0.2 

3 

232 

235 

0.38 

59.4 

40.5 

0.25 

3 

0.69 

45.1 

54.2 

0.80 

17 

0.85 

24.2 

74.3 

1.39 

19 
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TABLE  5.1  (Cont'd.) 


Temp. 

a 

*  c° 

b 

P. 

torr 

n=l 

-f- 

Mass  Spectrum  of  H  (HO) 

2  n 

2  3  4  5  6 

7 

Total  Ion 

Current 

_  14 

A  x  10 

1.55 

26.9 

71.2 

2.01 

9 

1.95 

24.6 

75.4 

• 

8 

2.85 

22.5 

77.5 

4 

232 

235 

3.50 

20.7 

79.3 

3 

335 

315 

0.85 

2.86 

89.6 

7.54 

90 

• 

2.20 

0.43 

84.9 

15.7 

130 

3.30 

0.82 

79.5 

19.7 

70 

5.43 

67. 8 

32.2 

15 

333 

315 

6.21 

65.1 

34.9 

11 

600 

493 

0.80 

78.6 

21.4 

55 

1.82 

59.5 

40.5 

50 

2.17 

56.2 

43.8 

30 

598 

512 

2.54 

52.3 

47.7 

20 

3.24 

46. 1 

53.9 

15 

595 

521 

4.23 

38.4 

61.6 

15 

* 

a 

is  the 

temperature 

of  the  thermocouple  in 

the 

ion 

source  bottom  while  b  is  the  temperature  of  the 
thermocouple  in  the  ion  source  top. 


. 


■ 
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Figure  5.1  Plots  of  log  K  versus  pressure  at  various 

X  /  z 


temperatures 
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P  (torn) 

Figure  5.2  Plots  of  log  K2  ^  versus  pressure  at  various  temperatures. 


t 


. 


119 


Figure  5.3  Plots  of  Log  K  versus  pressure  at  various 
temperatures . 


120 


Figure  5.4  Plots  of  log  ^  versus  pressure  at  various 
temperatures.  Plots  of  log  K  are  continued  on  Figure  5.5. 


■ 
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P  (torr) 


5.5  Continuation  of  lop  ^  plots  fiom  Fipuis  5.4. 
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Figure  5.6  Plots  of  log  K  versus  pressure  at  various 

Z)  f  O 


temperatures 
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temperatures  . 
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Figure  5.8 

Van't  Hoff-type  plots  of  the  equilibrium  constants.  0,  AMS;  0,  present  results. 
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5.2 


log  K 


n-l,n 


AH°  ,  +  AS° 

n-1  ,n  n-1 ,n 

2.3  RT  2.3  R 


AH0  1  at  298°K  is  obtained  from  the  slope  and  AS°  ,  at 

n  i,n  A  n-l,n 

298°K  is  obtained  from  the  intercept.  AG°  ,  at  298°K 

n- 1 ,  n 

results  from  extrapolation  of  the  straight  line  plots.  The 
values  for  these  thermodynamic  functions  are  listed  in 


Table  5.2. 


C.  Interpretation  of  the  Van't  Hoff  Plots 
Some  comments  are  in  order  on  the  reliability  of 
the  values  determined  by  the  van't  Hoff  plots.  First ,  because 

log  K  varied  rapidly  with  K,  the  error  in  the  slope  due  to 

1.  r  z 

the  uncertainty  in  the  temperature  measurement  (for  example, 
the  average  temperature  deviation  in  experiment  4  of  Table  5.1 
is  10°  and  in  experiment  5  is  50°)  leads  to  maximum  uncertainties 
AH0  ,  AG°  ,  and  AS0  of  ±15%,  ±20%,  and  ±25%  respectively. 
The  percentage  error  of  AH  was  calculated  by  drawing  error 
bars  for  the  temperature  on  the  van't  Hoff  plot  and  calculating 

AH°  from  the  steepest  and  shallowest  lines  which  could  be 

1,2 

drawn  through  the  error  bars  of  all  the  points.  The  errors 

in  AG°  were  obtained  by  extrapolation  of  the  lines 
298 

corresponding  to  the  maximum  and  minimum  values  of  AH°29g.  The 

errors  in  AS°  _  were  calculated  from  the  maximum  AH°  and 
298  298 

AG°298  and  minimum  AH°2gg  anc^  ^G°298  va-'-ues* 


in 
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TABLE  5.2 


Values 

for  the 

Thermodynamic 

Functions 

for 

the  Hydration  of 

the 

Pr< 

—Ah  ° 

-AC° 

-AS°  _  . 

n-1 ,  n 

i — 1 

i 

CO 

0^ 

CN 

1 

1 

n) 

298  ( 

n-1, 

n) 

(kcal/mole) 

(kcal/mole) 

(  e.u.) 

n-l,n 

Combined  AMS 

PMS* 

Combined 

AMS 

PMS* 

Combined 

AMS 

PMS 

0,1 

169** 

• 

- 

If  2 

36 

36 

25 

25 

33 

33 

2,3 

22.3 

22.3 

13.6 

13.6 

29 

29 

3,4 

17.2 

18.  3 

18.8 

8.5 

9.6 

9.5 

28 

32 

32 

4,5  . 

15.  3 

15.7 

15.3 

5.5 

5.6 

5.5 

33 

34 

33 

5,6 

13.0 

13.3- 

13.0 

3.9 

3.9 

3.9 

30 

31 

30 

6,7 

11.7 

11.7 

11.7 

2.8 

2.8 

2.8 

30 

30 

30 

7,8 

10.3 

10.3 

2.2 

2.2 

27 

27 

*  PMS  is  the  present  work, 


*  * 


From  Reference  (51  ) 


In  the  van't  Hoff  plot  for  K  the  line  has  been  drawn 
more  steeply  than  the  "best"  straight  line  through  the  points 
to  make  As°  fall  in  the  range  in  which  the  other  entropy 

values  fell.  AH°2g8,  AG02gg,  and  AS°2gg  for  the  "best"  straight 
line  are  numerically  smaller  by  1.7  kcal/mole,  0.8  kcal/mole, 
and  3  e.u.  respectively. 

Combining  the  results  from  the  two  mass  spectrometers 
has  both  an  inherent  advantage  and  disadvantage.  If  the 
deviation  of  the  results  from  the  best  straight  line 
through  all  of  the  points  is  due  to  random  measurement 
errors,  combining  the  points  from  the  two  instruments 
increases  the  accuracy  of  the  results.  However,  if  there 
is  a  systematic  error  between  the  two  sets  of  results, 
combining  the  points  increases  the  error.  In  the  K  plot, 
separate  lines  drawn  through  each  set  of  points  results 
in  a  AH0  value  of  -18.3  kcal/mole  for  the  AMS  and  -18.8 

0,4 

kcal/mole  for  the  PMS  (the  present  results) .  The  real 
value  could  be  between  -18.3  and  -18.8  or  it  could  be  -17.2 
kcal/mole  considering  the  line  through  the  combined  results. 
Fortunately  this  is  an  extreme  example  since  the  other 
results  in  which  both  instruments  were  used  to  measure 

K  .  have  better  agreement  either  separately  or  combined. 

n-1 ,  n 


.. 
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5 •  ^  Discussion  of  Results 

A.  The  Hydrates  as  a  Function  of  Pressure 

As  explained  in  section  4.3  the  relative  hydrate 
concentrations  can  be  calculated  by  use  of  equation  4.13  and 
4.14.  The  values  for  A  and  B  are  given  in  Table  5.3.  Figure 
5.9  illustrates  the  change  in  the  hydrate  concentrations  with 
pressure  at  temperatures  of  300°  and  400°.  This  type  of 
information  is  useful  to  radiation  chemists  studying  the  gas 

phase  radiolysis  of  a  polar  molecule.  For  example  Anderson 

« 

(  5  )  has  attempted  to  explain  a  temperature  effect  in  the 

y-radiolysis  of  water  vapor  by  the  variation  in  the  cluster 
size  with  temperature.  The  reactions  which  he  considers  are 

5.12  and  5.13.  The  variation  in  cluster  size  is,  therefore, 

5.12  H  0+  (H  O)  +  e“- - +H  +  (n+1)  Hn0 

3  2  n  2 

5.13  H30+  +  e~ - 2H  +  OH 

partly  responsible  for  changing  G (H^)  /  according  to  Anderson. 

B.  Entropy  Changes 

As  shown  in  section  4.3  a  value  can  be  theoretically 
calculated  for  the  entropy  loss  in  a  condensation  reaction. 
Comparison  of  the  experimental  and  theoretical  values  gives 
an  indication  of  the  validity  of  the  experimental  results. 
Since  the  calculations  depend  on  assumed  structures ,  it  is 
only  practical  to  calculate  the  simplest  case,  AS^  ^ 


The 


' 
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TABLE  5.3 


Constants 


for  the  Determination  of  K 
- - • — — — n-l,n 

Standard  State  1  torr 


n-1  ,n 


1/2 

2.3 

3.4 

4.5 


5,6 


6,7 


A 

B 

7.82 

10.1 

4.76 

9.05 

3.63 

8.86 

3.33 

10.0 

2.84 

9.46 

2.53 

9.24 

2.17 

8.50 

7,8 
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Figure  5.9 

Equilibrium  distributions  of  clusters  H  (H^O)^  Prec^cte<^  ky 
experimental  data;  logarithmic  scale  used  for  water  partial 
pressure  in  order  to  accomodate  a  wide  pressure  range. 

(a)  300 °K;  (b)  400°K.  Distributions  at  other  temperatures 

can  be  simply  calculated  from  data  in  Table  I. 
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translational  entropy  loss  by  the  Sackur-Tetrode  equation  is 
-33  e.u.  (See  section  4.3  for  an  example  of  this  type  of 
calculation) . 

Using  the  assumed  structures  for  H^O  and  H  (H^O)  ^ 
shown  in  Figure  5.10,  the  rotational  entropy  change  can  be 
calculated  from  the  type  of  equations  given  in  section  4.3. 
The  calculated  moments  of  inertia  are  also  given  in  Figure 
5.10.  The  resulting  theoretical  rotational  entropies  are 

-f*  *|* 

14.2  e.u.  for  HO  and  9.8  e.u.  for  H  (HO).  Using  the 
value  of  10.1  e.u.  for  H^O  calculated  in  section  4.3 


AS0  ^  .  =  -5.7  e.u. 

rot  1,2 

In  addition  to  the  overall  rotation  there  is  an 

internal  rotation  S°.  ^  ,  of  the  two  1F0  groups  in 

int  rot  2 

H+(H  0)  about  their  ion-dipole  axes.  The  corresponding 
entropy  term  can  be  calculated  from  the  partition  function 


for  free  rotation  (  48)  £>  by  equation  5.3  where  the  reduced 

Fr 


5.3 


-Fr 


=  1 
a 


'  8tt  3  I  k  T 
r 


moment  of  inertia  I  is  given  by  equation  5.4.  For  the 
complex  under  consideration  1^  equals  equals 


5.4 


I  =  Ii  Io 
r  12 


Il+I2 


-40 

the  moment  of  inertia  about  the  ion  dipole  axis,  1.91  x  10  g 


2 

cm  . 


4 


. 
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■I"  4” 

Figure  5.10  Assumed  structures  of  H^O  and  H  (H2°)2  for 
the  theoretical  calculation  of  AS0  . 
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On  substitution  in  equations  5.3  and  4.8,S°  =  +  4  7  e  u 

int  rot 

The  remaining  terms  needed  to  calculate  AS0 

i.  t  z 

are  vibrational.  The  magnitude  of  these  terms  is  expected 
to  be  very  small  because  of  the  large  experimental  values  of 
AH  and  AH  which  are  -169  (  51)  and  -36  kcal/mole 

f  J-  -L  /  Z 

respectively.  These  values  indicate  a  high  bond  energy.  From 
equation  4.11  the  torsional  vibrations  similar  to  those 

-f* 

described  for  K  (H^O)  are  expected  to  have  a  negligible  entropy. 

The  high  bond  energy  also  indicates  that  the  force  constants 

for  the  other  vibrations  are  high.  As  can  be  seen  from 

equations  4.10  and  4.11  this  results  in  a  small  value  for 

entropy  for  the  new  vibrations.  We  are  naively  assuming  the 

original  vibrations  to  be  unchanged.  Consequently  taking 

AS0  .  =0  qives  AS°  ^  =  -34  e.u.  The  difference  between 

vib  1/2 

the  theoretical  and  experimental  entropy  values,  +  1  e.u.y 
is  certainly  within  the  combined  errors  of  theory  and 
experiment . 

C.  Enthalpy  Changes 

The  values  of  the  thermodynamic  functions  for 
reaction  7,8  are  compared  to  the  values  for  the  condensation 
of  water  vapor  and  the  values  for  the  gas  phase  reaction  5.14 

5.14  H20  +  H20^  (H202 

in  Table  5.4.  The  fact  that  AH°?  Q  is  close  to  the  heat  of 
condensation  of  water  need  not  mean  that  the  effect  of  the 


■ 
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TABLE  5.4 


Comparison  of  Values  for  the  Thermodynamic  Functions 


corresponding  to  the  condensation  of  a  water  molecule 

onto  H+(H^O)._,  and  H^O  ^  and  into  a  pool  of  liquid 


water  with 

a  planar  surface. 

H+(H20)7(g) 

^  a 

H_0,  \ 

2  (g) 

_  b 

H2°(l) 

AH°298 

-  10.3 

-  5.0 

-  10.52 

kcal/mole 

AG°298 

-  2.2 

+  2.7 

-  2.05 

kcal/mole 

AS°298 

-  2.7 

-25.8 

-  28.39 

e.u. 

From  Reference  ( 59  )  • 


From  Reference  ( 61  )  • 
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ionic  change  is  negligible  beyond  the  octahydrate.  From 
the  enthalpy  given  in  Table  5.4  for  reaction  5.14,  it  is 
apparent  that  the  heat  of  condensation  for  a  water  molecule 
condensing  onto  a  small  electrically  neutral  gaseous 
nucleus  is  considerably  less  than  the  heat  of  condensation 
into  liquid  water  with  a  planar  surface. 

One  reason  for  this  is  that  the  enthalpy  of  condensation 
depends  on  the  radius  of  curvature  of  the  water  surface. 

Dufour  and  Defay  (62)  estimate  the  heat  of  evaporation 

O 

from  a  droplet  of  radius  10 A  to  be  0.94  of  that  from  a 

O 

plane  surface.  Extending  the  calculation  to  4A,  an  assumed 
value  for  the  radius  of  the  octahydrate,  one  obtains 

O 

AH  (r  =  4A) /AH  (r  =  °°)  =  0.85,  which  gives  8.8 

evap  evap 

kcal/mole  for  the  heat  of  evaporation  from  a  neutral  cluster 

O 

with  r  =  4A.  Since  the  calculation  depends  on  the  macro- 

O 

scopic  concept  of  surface  tension,  its  extension  to  a  4A 
radius  is  tenuous.  However,  the  calculation  does  indicate 
that  a  significant  fraction  of  AH7  q  is  made  up  of  ion- 
molecule  attractions  rather  than  molecule-molecule  attractions. 

D.  The  Structure  of  the  Complex  Ions 

Consideration  of  the  structure  of  H2n+1  0^*  involves 

two  models.  The  first  model  is  similar  to  that  proposed  for 

K+(h  0)  ,  namely  a  central  ion  with  coordinated  water  molecules. 

2  n 

The  central  ions  in  the  observed  H  x  C>n  ions  may  be  predominantly 

The  ion  with  molecular  weight  19  is 


hydronium  ions  or  hydrogen  ions. 


'  "  **  ¥***  n0  i+ n ^  sdi  n£  anoi  Xo^ireo  ®rfT 
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■{■  -j.  i 

undoubtedly  HO  rather  than  H  (H^O) .  The  notation  H  (H  0) 

^  2 

implies  that  all  the  ionic  charge  resides  on  one  hydrogen, 
which  is  incorrect  (63).  The  next  ion  in  the  series, 
molecular  weight  37,  can  be  represented  equally  well  as 

-f* 

H  (H^O)  2  °r  ^3^  (**2^  *  ^ne  Can  the  Pr°ton  as  being 

equally  bonded  to  the  two  water  molecules.  If  the  higher 

hydrates  are  of  the  form  HO  (HO)  ,  a  break  might  be 

3  2  to 

expected  in  the  thermodynamic  functions  between  reactions 
5.15  and  5.16.  The  reason  is  that  although  the  first  three 
water  molecules  can  coordinate  to  the  hydrogen  atoms  of  H^0+, 


5.15 

H  cT(H  O)  0 
3  2  2 

+  h2o 

-> 

-e 

H30+(H2°) 

5.16 

H  O*  (H  O)  _ 
3  2  3 

+  H2° 

-> 

-e 

h3o+(h2o) 

the  fourth  water  molecule  must  coordinate  in  a  new  way.  This 

way  may  be  the  formation  of  an  outer  complex  of  the  type  shown 

in  Figure  4.21  or  the  coordination  of  the  water  molecule  to  the 

oxygen  atom  of  the  hydronium  ion. 

Figures  5.11  and  5.12  show  AH°  and  AG°  at  298°K 

r  n-l,n  n-l,n 

as  a  function  of  n-l,n.  These  curves  do  not  show  any  obvious  breaks. 

For  this  reason  the  first  derivative  curves  of  AH°  .  and 

n-l,n 

AG°  are  also  plotted  in  Figures  5.11  and  5.12  respectively. 

n-l,n 

Although  A(AH°  ^  shows  a  slight  break  between  reactions  5.15 

and  5.16,  A  (AG°  .  )  does  not.  Although  the  possibility  that 

n-l,n 

H  0+  is  the  central  ion  is  not  ruled  out  since  there  is  a  slight 
3 

break  and  since  the  AH0.  _  value  from  the  present  work  is  somewhat 

4  #  5 

inaccurate  due  to  a  small  curvature  in  the  van't  Hoff  plot,  it  would 
that  H+  is  more  likely  the  central  ion. 


appear 
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Figure  5.11 

Plot  of  AH°  .  and  A  (AH°  .  )  vs.  n-l.n. 

n-l,n  .  n-l,n 

of  the  reactions  decrease  quite  regularly. 

break-  is  indicated  between  AH  and  AH,  _ 

3,4  4,5 


The  enthalpies 
A  possible 
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Figure  5.12 


Plot  of  AG°  .  and  A (AG°  _  )  at  298°K  vs.  n-l,n.  Standard 

n-l,n  n-l,n 

state  of  water  vapor  is  1  atm.  Decrease  of  free-energy  changes 


is  quite  regular  showing  that  stabilities  of  clusters  change 


quite  regularly. 
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This  suggests  a  second  model  without  a  fixed  central  ion. 
In  this  model  the  ionic  charge  moves  from  molecule  to  molecule 
in  the  cluster  by  proton  transfer.  This  model  has  as  its  basis 
the  high  proton  mobility  observed  in  liquid  water.  Such  a 

model  would  suggest  smooth  changes  in  AH0  ,  and  AG°  , 

^  n-l,n  n-l,n 

with  n-l,n  since  each  successive  condensation  of  a  water 

molecule  would  give  rise  to  similar  interactions. 

These  two  models  are  not  very  different  if  it  is 

considered  that  H+(H  0)  and  H_0+(H.0)  ..  are  important 

2  n  3  2  n-1 

transient  structures  in  the  mobile  charge  model. 


. 
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6  *  SOLVATION  OF  THE  AMMONIUM  ION  BY  AMMONIA  MOLECULES 

6. 1  Introduction 

The  study  of  the  NH.+(NH  )  system  was  undertaken  to 

4  3  n 

determine  AH0  ,  AG°  ,  ,  and  AS°  ,  for  reaction  6.1. 

n-l,n  n-l,n  n-l,n 


6.1  (n— l,n)  NH  +  (NH  )  .  •+  NH  NH  .+  (NH_) 

4  3  n-1  3  ■* —  4  3  n 

/ 

+ 

The  -system  NH  .  (NH_)  is  of  interest  because  of  its  similarity 
4  3  n 

to  H+(H^0)^.  In  Chapter  5  the  experimental  thermodynamic 

values  suggested  that  the  hydrated  species  being  investigated 

4"  4" 

was  H  (H„0)  rather  than  H_0  (H„0)  n.  Previous  work  by 
2  n  3  2  n-1 

Kebarle  (  15  )  on  reaction  6.1  for  steps  n-l,n  =  2,3;  3,4; 

and  4,5  indicated  the  solvated  species  was  NH.  (NH_)  .  The 

4  3  n 

structure  was  also  indicated  from  a  study  of  NH.4(H  0)  by 

4  2  n 

Hogg  and  Kebarle  (  17  ) .  The  present  work  was  done  to  obtain 
a  complete  set  of  thermodynamic  data  on  reaction  6.1  to 

-f" 

confirm  or  refute  the  structure  NH  (NH  )  . 

4  3  n 


6.2  Reactions  Leading  to  NH A  (NH0) 

■  ■  ■  - -  - - - — — ■  ■■ — — h — - j  n 

4_  -f- 

NH  and  NH  are  the  main  primary  ions  formed  m  the 
3  2 

proton  beam.  As  is  the  case  in  irradiated  water  vapor,  the 
primary  ions  react  rapidly  to  give  a  protonated  species. 


4~ 

+ 

6.2 

nh3 

+ 

nh3  - 

- NH 

4 

+  NH2 

+ 

6.3 

NH2 

+ 

nh3  - 

- ->  NH 

4 

+  NH 

The  thermal  energy  reaction  rates  of  reactions  6.2  and  6.3 


IS)  iBBTt&tif  arfr 
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are  5.2  x  10  and  8.4  x  10  cm  /molecule-sec.  respectively 
(  64  ) .  The  further  reactions  of  NH^  are  given  by  reactions 
6.4,  6.5,  and  6.1  where  M  is  a  third  body  to  remove  the 


energy  released  on  condensation.  The  equilibrium  constant 


6.4  (0,1)  NH4+  +  NH  JU  NH4+(NH3) 

+  M  + 

6.5  (1,2)  NH4  (NH3)  +  NH3 - NH4  (NH) 


K  is  given  by  equation  6.1  where  I  +,  ,  and  I  + 

n-l,n  NH,  (NHJ_  NH  ,  (NH,,) 


4  '  3  n 


4  '  3  n-1 


6.1 


K  =  INH  +(NH  ) 

n-l,n  4  3  n 


INH.+  (NH.)  ?NH 
4  3  n-/  3 


•f 

are  the  mass  analyzed  ion  currents  of  NH,  (NHJ  and  NH  „  (NH  ) 

4  3  n  4  3  n-1 


respectively.  P  is  the  pressure  of  ammonia  in  the  ion  source. 

NH3 


6.3  Effect  of  Ion  Exit  Slit  Dimensions  on  Results 

The  dimensions  of  the  ion  exit  slit  in  the  ion  source 
were  found  to  affect  the  Jesuits.  Initially  the  equilibrium 
constants  were  measured  at  temperatures  up  to  590 °C  with  a 
5y  x  4  mm  ion  exit  slit.  At  590°C  the  measurement  of  K 

U  f  1 

was  interfered  with  by  the  presence  of  an  impurity  which 
gave  an  ion  current  of  K+  which  was  seventy  times  larger 
than  the  total  ion  current  due  to  NH4+(NH3)^.  Consequently 
the  temperature  was  reduced  to  330 °C  and  the  constants 

K  and  K  were  measured.  These  values  were  erroneous. 

0,1  1,2 

K  was  a  factor  of  two  lower  and  K  ten  times  higher 

0  » 1  1/4 
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than  the  previous  results  predicted.  On  inspection 
of  the  ion  source  ion  exit  slit,  it  was  found  that  the 
slit  had  buckled  upon  cooling  from  590°C.  The  area  of 
the  buckled  slit  was  much  larger  than  the  original  area. 

A  new  ion  exit  slit  was  made  15y  x  1.7  mm.  The  equilibrium 
constants  measured  with  this  slit  installed  in  the  ion 
source  were  in  agreement  with  the  undistorted  5y  x  4  mm 
slit  as  shown  in  Figure  6.8. 

At  low  temperatures  the  results  using  the  15y  x  4  mm 
slit  became  unacceptable.  Figure  6.1  shows  the  constants 
measured  as  a  function  of  pressure  at  -32 °C.  At  a 
pressure  of  0.57  torr  an  increase  in  the  ion  acceleration 
potential  between  the  ion  source  and  the  first  acceleration 
plate  (draw-out  plate)  from  360  v  to  550  v  changed  K 

3  9  ^ 

from  6.0  to  1.4.  This  indicates  that  a  substantial 
proportion  of  the  mass  analyzed  ions  have  undergone  a 
collision  just  outside  the  ion  source.  The  higher  ion 
acceleration  causes  a  greater  number  of  ionic  clusters 
to  undergo  dissociative  collisions. 

To  reduce  the  number  of  collisions  outside  the  ion 
source,  the  ion  exit  slit  area  was  reduced  by  about  a 
factor  of  five.  The  new  dimensions  were  1.5y  x  3.5  mm. 

The  measurements  of  the  constants  with  this  slit  in  the 
ion  source  gave  pressure  independent  equilibrium  constants 
at  all  temperatures  as  illustrated  in  Figures  6.4,  6.5,  6.6, 
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P  (torr) 

Figure  6.1  Equilibrium  constants  affected  by  collisions 


outside  the  ion  source. 
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and  6.7.  The  constants  were  also  in  agreement  with  the 
previously  measured  constants  which  are  also  shown  in 
Figures  6.4,  6.5,  6.6,  and  6.7, 

6 . 4  Results 

Some  representative  mass  spectra  are  listed  in  Table  6.1 

Figures  6.2  -  6.7  Show  the  plots  of  K'  ,  versus  ammonia  pressure 

n-l,n 

obtained  at  temperatures  from  -57°  to  485°C.  Figure  6.8  shows  van't 
Hoff  type  plots  of  the  equilibrium  constant.  Values  for  the 

thermodynamic  functions  AH0  n  ,  AG°  n  ,  and  AS0  ,  at  298°K 

n-l,n  n-l,n  n-l,n 

obtained  from  equation  5.2  are  listed  in  Table  6.2.  Included  in 
Table  6.2  are  the  values  obtained  by  Hogg,  Haynes,  and  Kebarle  (  15  ) . 
The  free  energy  values  are  seen  to  be  in  good  agreement,  although 
the  enthalpy  and  entropy  values  are  not  in  good  agreement.  Table 
6.2  does  not  include  values  for  reaction  5,6  because  of 
insufficient  data. 

6 . 5  Discussion 

A.  NH  +  (NH.J  as  a  Function  of  Pressure 

- 4 — u—g—n— - - - — - - — ■ — — 

By  the  method  outlined  in  section  4.3  the  relative  ionic 
concentrations  can  be  calculated  by  use  of  equations  4.13  and  4.14. 
.The  values  for  A  and  B  are  listed  in  Table  6.3.  Figure  6.9 
illustrates  the  change  in  relative  concentrations  as  a  function  of 
ammonia  pressure  at  300°K.  As  explained  in  section  5.3  the  cluster 
concentrations  are  of  interest  to  radiation  chemists  studying  the 
vapor  phase  radiolysis  of  a  polar  molecule. 

B.  Entropy  Changes 

The  entropy  change  occurring  in  reaction  6.4  is 


"*  V 


. 
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Figure  6.2 


Plots  of  log  K 

U  /  1 


versus  pressure  at  various 


temperatures . 


0,  5  y  x  4  mm  ion  exit  slit  dimensions. 


. 


P  (t  o  r  r) 


Figure  6.3 
temperatures 


Plots  of 
0,  5  y 


L°g  K1j2 
x  4  mm; 


versus  pressure  at  various 
□  ,  15  y  x  1.7  mm  ion  exit 


slit  dimensions 
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P  (iorr) 

Figure  6.4  Plots  of  log  K?  3  versus  pressure  at  various 
temperatures .  0,  5  y  x  4  mm;  D , 15  y  x  1.7  mm;  /S  ,  1.5  y  x 


3 . 5  mm 


ion  exit  slit  dimensions. 


148 


P  (torr) 


Figure  6.5  Plots  of  Log  K  versus  pressure  at  various 

temperatures.  0,  5  y  x  4  nun;  □  ,  15  y  x  1.7  mm; 

/S  ,  1.5  y  X  3.5  mm  ion  exit  slit  dimensions. 
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Figure  6.6 
temperatures 


Plots 

A, 


of  log  K  versus  pressure  at  various 
1.5  y  x  3.5  mm  ion  exit  slit  dimensions. 


Log  K- 
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Figure  6.7 
temperatures . 


Plots  of  log  K_  versus  pressure  at  various 

J  /  6 

A,  1.5  y  x  3.5  mm  ion  exit  slit  dimensions. 
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Figure  6.8  Van ' t  Hoff- type  plots  of  the  equilibrium 
constants.  0,  5  y  x  4  ram;  D  ,  15  y  x  1.7  ram;  A,  1.5  p  x 
1.7  mm  ion  exit  slit  dimensions. 


100 
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TABLE  6.1 


Representative 

Data 

Given  as  the 

Mass 

Spectra 

in 

Per 

Cent 

Total  Ionization 

Temp . 

*  c° 

P. 

Mass  Spectrum  of  NH^+ 

(NH  ) 

3  n 

Total  Ion 

Current 

14 

A  x  10 

No. 

a 

b 

torr 

n=0 

1 

2 

3 

4 

5 

1 

-2 

-5 

0.30 

1.1 

50 

49 

1 

0.63 

0.8 

30 

69 

3 

0.79 

1.3 

26 

73 

0.21 

4 

-2 

-6 

0.98 

1.6 

24 

74 

0.45 

4 

1.22 

2.0 

27 

70 

0.60 

3 

-2 

-6 

1.85 

2.3 

27 

70 

1.0 

1 

2 

93 

105 

0.20 

24 

71 

4.4 

1 

0.40 

8. 1 

81 

11 

4 

1.00 

4.4 

74 

21 

0.07 

10 

1.41 

.7 

68 

31 

0.14 

10 

2.22 

1.0 

59 

40 

0.34 

9 

88 

100 

3.21 

2.1 

52 

45 

0.61 

5 

3 

379 

391 

0.20 

95 

5.0 

0.  3 

0.67 

74 

26 

0.027 

2 

1.  31 

59 

41 

0.076 

10 

379 

391 

2.23 

46 

54 

0.16 

40 

3.04 

40 

60 

0.27 

100 

381 

398 

4.36 

36 

64 

0.30 

100 

*a  is 

the 

temperature 

of  the  thermocouple  in 

the  bottom 

of  the  ion  source  while  b  is  the  temperature  of  the 
thermocouple  in  the  ion  source  top. 


. 
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TABLE  6.2 


Experimental  Thermodynamic  Values  for  the  Gas  Phase  Reactions: 


NH 

4 

+  (NH  )  + 

3  n-1 

NH - > 

3  .< - 

NH  /(NH  ) 

4  3  n 

i — ( 

i 

at 

298°K  and 

1  atm. 

n-1  ,n 

-AH°  , 
n-1  ,n 

-AG°  , 
n-1  ,n 

-AS°  , 
n-l,n 

Present 

(kcal/mole) 

Present 

(kcal/mole) 

Present 

(e .  u. ) 

work 

Ref.  15 

work 

Ref.  15 

work. 

Ref.  ; 

0,1 

27 

17.5 

32 

1,2 

17 

9.0 

27 

2,3- 

14.0 

17.8 

6.1 

6.4 

27 

38 

3,4 

13.1 

15.9 

3.6 

3.8 

32 

40 

4,5 

7.5 

(9)b 

0.2 

0.5 

25 

(33)  a 

a  Estimated. 

b  Obtained  from  AG°  and  the  estimated  AS° 

hi  /  -> 


V 
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TABLE  6.3 


Constants  for  the  Determination  of  K 
Standard  state  1  torr. 


n-1  ,n 


Reaction  n-l,n 

A 

B 

0,1 

5.87 

9.84 

1,2 

3.69 

8.75 

2,3 

3.05 

8.75 

3,4 

2.85 

9.84 

4,5 

1.63 

8.  31 
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expected  to  be  similar  to  the  corresponding  reaction  in  the 
H  (H^O)^  sys^em/  reaction  5.3.  The  values  for  respective 
experimental  entropy  changes  are,  in  fact,  very  close,  -32 
and  -  33  e.u.  This  agreement  plus  the  theoretical 
calculation  of  the  entropy  change  for  reaction  5.3  suggests 


the  experimental  entropy  value  for  AS°^  ^  for  reaction  6.4 
is  approximately  correct. 


C.  Comparison  of  Results  to  Other 
Ammonia  Condensation  Reactions 

The  values  for  the  thermodynamic  functions  of  reaction 
4,5  are  compared  to  the  values  for  reactions  6.6  and  6.7  in 


6.6 


NH  ,  +  NH  .  T 
3(g)  3(g) 


*  (NH3>2(g) 


NH 


6.7 


NH 


3U) 


3(g) 


*  NH 


3U) 


Table  6.4.  The  magnitude  of  the  entropy  changes  are  seen  to 
be  similar  for  all  three  reactions.  Because  of  the  low  heat 
of  condensation  of  ammonia,  the  forces  in  the  6th  cluster  are 
predominantly  ion-molecule  rather  than  molecule-molecule 


forces . 


D.  The  Structure  of  the  Solvated  Ions 

As  discussed  in  the  previous  chapter,  there  is  a  choice 

between  two  models  for  the  structure  of  the  solvated  ions.  In 

.  + 

the  first  model  a  central  ion,  either  NH^  or  H 


/ 


is  surrounded 


' 
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TABLE  6.4 

Comparison  of  Values  for  the  Thermodynamic  Functions 

of  Condensation  Type  Reactions 
at  298°K 


Reaction 

AH0 (kcal/mole) 

AG° (kcal/mole) 

AS  °  (e.u 

4,5 

-7.5 

-0.2 

-25 

6.6a 

-4.4 

+3.5 

-26.8 

6.7b 

-4.75 

+2.17 

-23.2 

a.  Ref.  (59)  . 


b.  Ref.  (61)  . 
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+ 

by  NH^  molecules.  If  the  central  ion  is  NH^  ,  one  might 
expect  a  break  in  the  thermodynamic  functions  upon  addition 
of  the  fifth  ammonia  molecule.  The  reason  is  that  only 

four  ligands  can  be  coordinated  with  the  hydrogen  atoms  of 

+ 

NH4  .  Consequently,  the  addition  of  the  fifth  ammonia  molecule 

requires  either  the  beginning  of  an  outer  shell  of  ammonia 

molecules  or  a  rearrangement  of  the  complex,  possibly  to 

H+ (NH  )  .  As  can  be  seen  from  Table  6.2,  the  values  for 

3  5 

all  three  thermodynamic  functions  take  a  sharp  drop  in 

absolute  magnitude  between  reactions  3*4  and  4,5.  This 

supports  the  idea  that  the  central  ion  of  the  complex  for 

n=l  to  4  is  NH4+.  An  analysis  of  the  relative  potential 
+  + 

energies  of  H  (NH  )  and  NH  (NH  )  (NH  )  would  be  required 

3  5  4  3  4  3 

to  determine  the  central  ion. 

We  may  now  consider  why  the  ammonia  system  has  the 
central  ion  NH  +  while  the  water  system  (section  5.3  D) , 
apparently  has  either  the  central  ion  H  or  an  ion  H 
which  is  free  to  move  about  in  the  cluster.  The  central 
ion  H+  possibility  for  the  ammonia  cluster  cannot  be 
ruled  out  on  theoretical  grounds  without  a  detailed  study. 
However,  the  mobile  H+  possibility  is  less  likely  in  the 
ammonia  system  than  in  the  water  system  since  the 

"h 

phenomenon  of  the  anomalously  high  mobility  of  H  which 
occurs  in  liquid  water  has  been  found  absent  in  liquid 
ammonia  in  an  investigation  by  Gill  (  65  )  . 
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7 .  CONCLUSIONS  FROM  THE  GAS  PHASE  SOLVATION  RESULTS 

7 . 1  Introduction 

The  reason  for  performing  this  work  was  to  obtain 
information  on  gas  phase  solvation  reactions.  It  was 
hoped  that  the  gas  phase  solvation  results  would  be 
useful  to  researchers  studying  solvation  phenomena  in  the 
liquid  phase.  The  advantage  of  the  present  work  over 
liquid  phase  work  was  that  the  empirical  formula  and 
relative  concentration  of  each  solvated  species  could  be 
determined  directly  by  mass  analysis. 

7 . 2  Structure  of  the  Solvated  Ions 

The  observation  of  solvated  species  in  the  present  work 
proved  the  existence  and  stability  of  certain  ionic  clusters. 
From  the  smooth  change  in  the  thermodynamic  values  for  AH0  , 

AG°  ,  and  AS°  n  for  the  reaction  7.1,  it  was  concluded 

n-l,n  n-l,n 

7 . i  H+(H20)n-l  +  H2°  < *  H+(H20)n 

in  section  5.3  D  that  notation  H  (H^O)^  represents  the 

ionic  clusters  than  H30+  (H20)  n_r  The  basis  for  conclusion 

is  the  assumption  that  if.  the  cluster  had  a  central  ion  H.^0  ,  a 
break  would  occur  in  the  thermodynamic  values  between  reactions 

3,4  and  4,5. 

Similarly ,  a  central  ion  NH^  would  be  indicated  for  the 

cluster  NH  + (NH  )  if  a  break  occurred  in  the  thermodynamic 
4  3  n 


. 
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values  between  the  specific  reactions  3,4  and  4,5  of  the  general 
reaction  7.2.  In  Chapter  6  a  definite  break  was  found  at  this  point. 


7.2  NH  h(NH  )  +  NH  - *  NH  +(NH  ) 

4  3  n-1  3  -  4  3  n 

The  question  of  the  identity  of  the  central  ion  does  not 

arise  in  the  hydration  of  the  potassium  ion.  Consequently  it  might 

be  expected  that  the  plots  for  the  thermodynamic  values  for  reaction 

7.3  as  a  function  of  n-l,n  would  be  completely  structureless.  However, 


7.3 


K+(H  0)  +  H  O - »  K+(H  0) 

2  n- 1  2  -**-  2  n 


as  seen  in  Figure  4.20  this  is  not  the  case.  As  indicated  in  section 

4.4  the  "bump"  in  the  AH°  ,  plot  versus  n-l,n  may  be  due  to  the 

n-l,n 

non-uniform  changes  in  the  ligand- ligand  distances  with  increasing 
cluster  size.  This  affects  the  changes  in  the  dipole-dipole 


repulsion  energies  which  are  primarily  responsible  for  the  change 

O 


of  AH  ,  with  n-l,n. 
n-l,n 

We  may  now  consider  why  the  plots  of  the  thermodynamic  values 

for  reaction  7.3  show  some  structure  while  the  thermodynamic  values 

for  reaction  7.1  produce  a  smooth  curve.  The  relative  lack  of 

structure  for  the  thermodynamic  values  for  reaction  7.1  is  an 

indication  that  the  proton  and  its  surrounding  water  molecules  are  not 

+ 

held  in  fixed  positions  while  K  is  most  probably  fixed.  Instead 
the  charge  in  H+(H20)^  can  move  about  the  cluster  by  proton 
transfer. 

A  second  means  of  determining  the  structure  of  the  ions 
is  from  a  plot  of  the  weighted  average  number  of  ligands 


::<w  3c  Jolq  4  noil  ei 
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per  ion  as  a  function  of  pressure.  Such  plots  are  illustrated 

in  Figure  7.1A  and  7.1B.  For  K+(H  0)  and  H+ (H  0)  no 

2  n  2  n 

single  ionic  species  appears  to  exhibit  dominant  stability 
in  the  pressure  range  considered.  However,  the  plot  of 

-f*  #  i 

NH4  (NH3)n  indicates  a  high  relative  stability  of  NH^  (NH3) 

-3  -1 

over  the  pressure  range  10  to  10  torr. 

One  conclusion  from  the  above  considerations  is  that  if 
a  single  structure  such  as  as  proposed  by  Eigen  and 

de  Maeyer  (66  )  has  special  stability  in  liquid  water,  the 
special  stability  is  created  by  the  effect  of  the  bulk 
solvent  on  the  ionic  cluster. 

7.3  Effect  of  the  Liquid  Phase  on  the  Ionic  Clusters 

A.  Ligand  Bond  Energies  in  the  Liquid  Phase 

One  way  in  which  the  present  AH°  ,  values  can  be 

n-l,n 

related  to  the  condensed  phase  is  through  the  theoretical 

arguments  of  Bernal  and  Fowler  (  67  )  .  They  have  proposed  the 

following  relating  to  the  ionic  solvation  in  liquid  water. 

The  hydration  of  an  ion  will  occur  when  the  potential 

energy  of  a  water  molecule  coordinated  to  an  ion  is  less  than 

that  of  a  free  water  molecule.  The  free  molecule  has  four 

neighbors  while  the  ion-bonded  molecule  has  three,  an  ion  and 

two  water  molecules.  The  condition  for  an  ion  having  at 

least  one  coordinated  water  is,  therefore,  that  the  potential 

energy  AE°  for  a  water  molecule  due  to  coordination  to 
n(£) 
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P  (torn) 

Figure  7 . 1A  Plot  of  the  weighted  average  number  of  ligands  per  ion  versus  pressure. 
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Log  P 

Figure  7. IB  Plot  of  weighted  average  number  of 
ligands  per  icn  versus  log  pressure.  The  log 
scale  is  used  to  cover  a  wide  pressure  range. 
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an  ion  be  less  than  the  potential  energy  2  P.E.  for  a 

H .  B . 

water  molecule  due  to  coordination  to  two  water  molecules. 

2  P.E.  is  calculated  by  Bernal  and  Fowler  to  be  -15 

n  .  d  . 

kcal/mole.  This  fact  gives  rise  to  equation  7.1  for  the 
ion-water  molecule  bond  energy  B.E. 

7.1  B.E.  =  AE°  ...  +15 

n  (£) 

AE°  ...  can  be  set  equal  to  AE°  ,  , . the  potential  energy 
n(£)  ^  n(g)'  * 

for  a  water  molecule  due  to  coordination  to  an  ion  in  the  gas 
phase.  This  is  not  unreasonable  since  Bernal  and  Fowler  appear 
to  have  taken  the  dielectric  constant  as  unity  in  their 
calculation  of  P.E. 

H  •  B  • 

Since  in  section  4.4  it  was  found  that  AH°  .  y  AE°  .  .  , 

n-l,n  ^  n(g) 

O 

AH  ^  AE°  , „ , .  Equation  7.1  can  be.  rewritten  as  equation  7.2. 

n-l,n  ^  n(£)  H 

7.2  B.E.  =  ’ AH°  .  +  15 

n-l,n 

The  bond  energies  calculated  for  K  (^0)^  —  ^2^  &n<^  ^  ^H2^^n  —  ^2^ 

with  the  AH°  values  from  Table  4.2  and  5.2  respectively  are 

n-l,n 

■f 

listed  in  Table  7.1.  The  bond  energies  for  K  (H2°)n  —  H2°  show 
that  the  potassium  ion  in  solution  may  have  one  or  two  coordinated 
water  molecules. 

However,  the  potassium  ion  is  considered  to  be  unhydrated 
in  solution  (55,68).  The  discrepancy  may  well  be  due  to 
approximations  in  the  theoretical  calculations.  Nevertheless, 
the  present  results  suggest  rather  weak  bonds  for  the  liquid 


' 


.  ,  . 
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TABLE  7.1 


I on- Water  Molecule  Bond  Energies  in  the  Liquid  Phase 


Bond 


K+  -  H20 
K+(H20)  —  H20 

K+(H20)2  “H2° 


Bond  Energy 
(kcal/mole) 

-  2.9 

-  1.1 
+  1.8 


h+(h2o)  — h2o 
h+(h2o)2-h2o 

H+ (H20) 3  H2° 

H+(H20)4  —  H2° 
H+(H20)5~H2° 


-  11 

-  7 

-  2 
0 


+  2 


% 

«• 
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phase  coordination  of  a  water  molecule  to  the  potassium  ion. 


The  bond  energies  for  H+(H  0) 

2  n 


H^O  indicate  that  the 


largest  expected  ionic  cluster  in  liquid  water  is  H  (H  0)_. 

2  5 

Because  of  the  relatively  strong  bond  energy  of  the  second 

molecule  and  the  zero  bond  energy  for  the  fifth  molecule, 

+  + 

the  dominant  ion  in  solution  could  be  H  (H  O)  or  H  (H  O)  . 

^  *3  b  A 


B.  Comparison  of  Experimental  Gas  Phase 

Free  Energy  Values  with  the  Free  Energy  of  Hydration 

The  free  energy  of  hydration  of  an  ion,  AG°  is 

hyd 

defined  as  the  free  energy  change  occurring  in  reaction  7.4  a  or  7.4  b 


7.4  a  M+.  ,  +  °°  HO...  — »  M+  .  +  00  H  O 

(g)  2  (£)  (aq.)  2  (£) 


7.4  b  X  .  +  00  H  0  - 

(g)  2  (£) 


X  +  w  H  0 

(aq.)  2  (£) 


The  determination  of  AG°,  ,  has  been  the  subject  of  a  number 

hyd 

of  papers.  The  present  gas  phase  free  energies  of  hydration 

AG°  ,  can  be  directly  compared  with  AG°,  ,.  In  a  recent 
n-l,n  hyd 

paper,  Kebarle,  Arshadi,  and  Scarborough  (21  )  have  compared 

the  values  of  AG°3  4  (measured  in  this  laboratory)  to  AG°h^d 

for  the  ions,  Na+,  K+,  F~,  Cl",  Br  ,  and  I  .  Although  this 

discussion  is  mostly  a  review  of  the  Kebarle  &~b.ciZ.  paper, 

the  present  approach  is  somewhat  different. 

If  one  assumes  that  both  AG°3  4  and  AG°hyd  have  the 

same  dependence  on  the  ionic  radius,  the  relative  values  of 

AG°  and  the  relative  values  of  AG°  should  be  similar. 

3,4  nYa 


. 
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On  this  basis  a  comparison  between  relative  free  energy  values 

is  given  in  Table  7.2.  The  values  for  AG°  in  column  B 

hyd 

resulted  from  Volta  potential  difference  measurements  between 
phases  separated  by  gas  (  69  ) •  The  values  in  columns  C  to  E 
result  from  the  Born  equation  (  70  )  ,  equation  7.3  in  which  N  is 

Avogadro's  number,  ze  is  the  ionic  charge,  r  is  the  ionic 
radius  and  D  is  the  dielectric  constant. 

7-3  -  AG°hya  “  f1  -  I 

2r  \  D 

The  experimental  radii  of  Gourary  and  Adrian  (71  ) 
result  in  the  best  theoretical  agreement  with  the  gas  phase 
free  energies.  The  values  of  Randles  (69  )  are  in  good 
agreement  with  both  those  values  from  Gourary  and  Adrian's 
radii  and  the  gas  phase  free  energies. 

7.4  A  Suggestion  for  Further  Research 

Because  the  gas  phase  solvation  reactions  are  relevant 
to  liquid  phase  solvation,  radiation  chemistry  (  6  ) ,  and  upper 
atmospheric  equilibria  (  4  ) ,  it  is  proposed  that  further  work 
be  carried  out  in  this  area.  One  interesting  project  would  be 
the  study  of  the  effect  of  the  ionic  charge  on  solvation.  Two 
reactions  which  might  be  studied  are  7.5  and  7.6.  The  reason 

7.5  Hg+(H20)n-l  +  H2°  Hg+(H20)n 

Hg+2(H20)n-l  +  H2°  ^  Hg+2(H20)n 


7.6 


d.V 
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TABLE  7.2 

Relative  Free  Energies  of  Hydration 


Ion 

AG03,4/AG°3,4(K+) 

A 

B 

AG°  VAG° 
hyd 

C 

hyd (K+) 

D 

E 

+ 

Na 

1.30 

1.22 

1.22 

1.27 

1.41 

+ 

K 

1.00 

1.00 

1.00 

1.00 

1.00 

F~ 

1.38 

1.23 

1.55 

1.28 

0.98 

Cl" 

0.91 

0.88 

1.14 

0.91 

0.74 

Br 

0.74 

0.81 

1.06 

0.83 

0.65 

i" 

0.51 

0.71 

0.95 

0.73 

0.59 

A.  From  data  compiled  in  Reference  (21)  . 

B.  From  Reference  (69)  experimental. 

C,D,E.  From  Born  equation  with  radii  from  References  (72,71,45) 


respectively . 
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+ 


+2 


for  proposing  mercury  ion  solvation  is  that  both  Hg  and  Hg 
can  be  created  by  ionization  of  Hg  vapor  by  the  100  kev 
proton  beam  which  is  available  in  this  laboratory.  One 
difficulty  in  this  study  is  that  reaction  7.7  may  interfere 
with  reaction  7.6. 


7.7 


Hg 


+2 


AH  =  -  6 .09  e.  v.  (  51  ) 
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1. 


2. 


3. 


4. 


5. 


6.  a 


b 


7. 

8. 


9. 


10 ,  a 


b 
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